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We review the fermionic Chern-Simons field theory for the Fractional Quantum 
Hall Effect (FQHE). We show that in this field theoretic approach to the problem 
of interacting electrons moving in a plane in the presence of an external magnetic 
field, the FQHE states appear naturally as the semiclassical states of the theory. 
In this framework, the FQHE states are the unique ground states of a system of 
electrons on a fixed geometry. The excitation spectrum is fully gapped and these 
states can be viewed as infrared stable fixed points of the system. It is shown that 
the long distance, low energy properties of the system are described exactly by 
this theory. It is further shown that, in this limit, the actual ground state wave 
function has the Laughlin form. We also discuss the application of this theory to 
the problem of the FQHE in bilayers and in unpolarized single layer systems. 



1 Introduction 

The fermionic Chcrn-Simons field theory for the FQHE0 was motivated by- 
Jain's observation that the Laughlin wave function could be reinterpreted as 
the wave function of charge-flux composites filling up a lowest Landau level, 
(i. e., a Vandermonde determinant), with each unit of charge being bound 
to an even number of flux quanta. Following Jain_.we will refer to these 
charge-flux composite particles, as composite fermionsua In Jain's picture, the 
electrons "nucleate" flux to screen enough of the external magnetic field, so 
that the "composite fermions" exactly fill an integer number of Landau levels 
associated with the unscreened part of the field. In this formulation, the FQHE 
is an Integer Quantum Hall Effect (IQHE) of the bound states. 

The starting point of Jain's approach was the Laughlin wave function. In 
1989, he proposed to write it in the suggestive factorized form 



*i(zi,...,ziv) =IJ( Zi ~ z ^ m 1 Xi{zi,---,Zn) (1) 



i<j 
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where xi is the wave function for a completely filled lowest Landau level 



N N 2 

Xi(z!,...,z N ) = Yl(zi-Zj) e *P(-J2^p)- (2) 

The phases associated with the first factor in Eq. |l| represent an even number 
(to — 1) of fluxes that are attached to each coordinate Zj where an electron is 
present. It is a crucial feature of this picture that the electrons bind to an even 
number of flux quanta and, in this way, they retain their fermion character. 
This observation gave rise to the picture of the FQHE as a ground state of 
"composite fermions" , where a composite fermion is an electron bound to an 
even number of fluxes. In this formulation, the main role of the long range 
correlations is to make it possible for the electrons to "nucleate" flux, and the 
FQHELcan be interpreted as an Integer Quantum Hall Effect of the bound 
states ii. 

The question that naturally arises here is how does a system of electrons in 
an external magnetic field manage to turn the correlations, which result from 
the electron-electron interactions, into "nucleated fluxes". In this Chapter 
we show that the key to. the answer is the Chern-Simons gauge theory. In 
1982 Wilczek's observeda that a particle current coupled to a Chern-Simons 
(CS) gauge field produced states with fractional statistics through the binding 
of particles to fluxes. Thus, if we are to get the Laughlin wave function by 
attaching (m — 1) fluxes to each electron, as suggested by Jain, it is natural 
to guess that the "right" theory must contain fermions (electrons) coupled 
to Chern-Simons gauge fields with an appropriate value of the Chern-Simons 
coupling constant 9. What is less clear is the origin of such a Chern-Simons 
gauge field for the problem of interacting electrons in a magnetic field. -In 
Quantum Field Theory, where the CS gauge theory was first introduced BEI, 
the Chern-Simons term in the action originates from the Parity Anomaly of 
relativistic fermions in 2+1 dimensions. Clearly, the electrons which live in 
quasi two-dimensional electron gases, are not relativistic. Also the real three 
dimensional system one is dealing with, neither breaks Parity by itself nor as a 
result of the presence of the device. The Chern-Simons gauge field must then 
be a result of the dynamics of interacting two dimensional electrons in the 
presence of an external magnetic field. For our purposes, the most important 
feature of the Chern-Simons term is not its relativistic invariance, but the fact 
that it is the only local gauge invariant theory which yields bound states of 
particles and fluxes. 

In this Chapter we derive a field theory for the FQHE in the fermion lan- 
guage in which the Chern-Simons gauge field appears explicitly in the problem 
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of interacting electrons in a magnetic field. We do so by first considering a 
theory in which the electrons, in addition to their mutual interaction, are cou- 
pled to both an external electromagnetic field and a Chern-Simons gauge field. 
In doing so, we use the fact that if the coefficient of the Chern-Simons term 
is chosen in such a way that an even number of flux quanta get attached to 
each electron, all the physical amplitudes calculated in this theory are identical 
to the amplitudes calculated in the standard theory, in which the CS field is 
absent (see reference El) . We further show that the Laughlin state is the semi- 
classical approximation of this theory. In the "classical" (i e., mean- field) 
approximation we get a picture very close to the one proposed by Jain. This 
theory, not only explains where do the fluxes come from, but also allows for 
the systematic calculation of corrections around this state in a manner similar 
to the semiclassical approximation in quantum mechanics and to the Random 
Phase Approximation (RPA) in many-body theory. A conceptually important 
feature of the Fermionic Chern-Simons theory is that its microscopic fields are 
composite fermions. Hence, at the level of mean field theory, it describes an In- 
teger Quantum Hall (IQH) state of composite fermions. Nevertheless, quantum 
fluctuations about this mean field state change this picture in a qualitatively 
fundamental way. They not only give the correct FQHE but also change the 
charge and statistics of the quasiparticles which become fractionally charged 
anyons, in agreement with the predictions of Laughlin's theory. The reason 
behind this behavior is the fact that the Chern-Simons fields require a local 
relation between charge and flux. This relation is replaced at the level of mean 
field theory by an average. The fluctuations correct this effect, forcing charge 
fluctuations to be accompanied by the necessary flux fluctuations. At long 
wavelengths, the leading gaussian fluctuations enforce this flux-charge rela- 
tion exactly, leading to the exact values of the Hall conductance, charge and 
statistics of the excitations as well as to the exact saturation of the sum rules. 
However, at shorter distances and higher energies, it is necessary to go beyond 
gaussian fluctuations to account for the correct correlations. 

This Chapter is organized as follows. In Section ^| we derive the Fermionic 
Chern-Simons theory of the FQHE for a single layer fully polarized two di- 
mensional electron gas (2DEG). In Section || we discuss the electromagnetic 
response functions and the spectrum of collective modes. In Section ^ we 
generalize this theory to the problem of the FQHE in bilayers and partially 
polarized and spin singlet FQH states in single layers. In Section [j| we use 
the Fermion Chern-Simons theory to derive the asymptotic long distance form 
of the wave functions of the FQH states. Finally, in Section |6] we present a 
summary of our results. 
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2 The Fermionic Chern-Simons Theory for the Fractional Quantum 
Hall Effect 

Consider a system of N electrons moving on a plane in the presence of an 
external uniform magnetic field B perpendicular to the plane. The electrons 
will be assumed to have an interparticle interaction governed by a pair potential 
y(|r|), for two electrons separated a distance \f\ on the plane. In this section 
we will assume that the magnetic field is sufficiently large so that the system is 
completely polarized and thus we will can ignore the spin degrees of freedom. In 
Section ^ we will consider polarization effects. The eigenstates ^(xi, . . . , xn) 
are eigenfunctions of the (first quantized) Hamiltonian H 

N 1 2 

6 = £<2m - - c m^)) + e Mm + £ - ( 3 ) 

i— 1 i<j 

where we have included the coupling to both the electromagnetic vector poten- 
tial A and scalar potential A$. Hence, we are dealing with N spinless fermions 
of charge — e and mass M. 

Our goal is to show that this system is equivalent to the same system 
but coupled to an additional statistical vector potential (fj, = 0, 1, 2) whose 
dynamics is determined by the Chern-Simons action, Sqs 

S cs = J d 3 x ie^T^. (4) 

for a suitably chosen value of 9. In Eq. || xq, x% and X2 represent the time and 
the space coordinates of the electrons respectively, and T vX is the field tensor 
for the statistical gauge field 

T vX = d v a x - d x a v . (5) 

The equivalent theory has a Hamiltonian H' which is identical to H , given in 
Eq. ||, except for the fact that the electrons are also coupled to the statistical 
vector potential a^. This is accomplished simply by setting 

-A^-A + a: eA -» eA + a . (6) 

c c 

The Chern-Simons gauge theory is a topological field theory. This means that 
the expectation values of the observables in this theory do not depend on the 
distance (i. e. , the metric) between the points of space-time on which they 
act, but only on topological properties. Being a gauge theory, its observables 
are path ordered exponentials of the circulation of gauge fields on closed loops 
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(Wilson loops! an their expectation values are functions of the linking numbers 
of these loopsQ. Since it is a topological theory, the action of the Chern-Simons 
gauge fields alone does not involve any energy scales. In fact, it is easy to see 
that the Hamiltonian is actually equal to zero (in the abscence of coupling to 
matter) . We will now follow standard arguments^ and show that the content of 
the Chern-Simons action reduces to a constraint and to a set of commutation 
relations for the gauge fields. Indeed, let us expand the Chern-Simons action, 
coupled to an external matter current jfj,{x), in components to read 

S cs = J d 3 x {[-j (x) + 9B(x)]ao(x) + 6a 2 (x)d a 1 (x)+a(x)-j(x)} (7) 

where j (x) is the particle density, j(x) is the matter current and B is the 
statistical flux. By direct inspection of Eq. |?], we see that the time component of 
the field ao (x) acts like a Lagrange multiplier field which imposes the constraint 

jo(x-) = 6B(x). (8) 

which is the analog of Gauss' law for this theory. At the quantum level, Eq. || is 
an operator constraint which selects the physical space of states i. e. , the gauge 
invariant states are made of charge-flux composites. The second term of the 
Chern-Simons action Eq. ^ implies that the momentum canonically conjugate 
to the field ai is just Qa-i and that, quantum mechanically they have to obey 
the equal-time commutation relations 

[a 1 (x),a 2 (x')]= t -5( 2 \x-x'). (9) 

u 

Finally, the last term in the action of Eq. [?], says that the Hamiltonian is just 

H ' = - J d 2 x j(x) ■ a(x) (10) 

which vanishes in the abscence of matter currents. 

Thus, for arbitrary values of the Chern-Simons coupling constant 0, the 
physical states are charge-flux composites: every particle of charge 1 carries 
a gauge flux equal to 1/9. Intuitively we expect that the wave functions for 
these composite particles should exhibit an Aharonov-Bohm effect which can 
be regarded as a change of statisticsQtj. In Chern-Simons theory this statistical 
transmutation is seen in the following way. Imagine that we have N particles 
which are sufficiently heavy so that their Feynman path integrals are dominated 
by reasonably smooth, well separated paths in 2 + 1-dimensional space-time. 
Consider for simplicity the case of two particles and let Ti and 1^ be their 
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respective paths. Let us denote by T = TilJ^ their joint paths. Let us 
further define a current j^{x) which carries one unit of charge and with support 
on r. A typical amplitude for this 2-particle system is given by the Wilson 
loop expectation value ErEl 

where the expectation values are calculated in the Chern-Simons theory and 
i^[r] is the linking number of the path T. From here it follows that pro- 
cesses involving the exchange of one pair of particles differ by the phase factor 
exp(i<5Ai^[r]), where we have defined the statistical angle S — ^ and Az^[T] 
is the change of linking number. For an exchange process, Ai/i[r] = ±1 and 
the phase factor is just exp(±i<5). Thus, S is regarded as the change (or shift) 
of statistics. By generalizing these arguments for an A-particle system we 
conclude that a system of fermions coupled to a Chern-Simons gauge field 
with coupling constant 9 behaves like a system of am/ons with statistical angle 
5 = 5* , measured with respect to Fermi statistics Eil. If 9 — J- ^- , where s is 
an arbitrary integer, then 5 = 2its and the system still represents fermions. 

In reference El we presented a detailed proof of the physical equivalence of 
two theories of particles coupled to a Chern-Simons gauge field with coupling 
constants 9 and 9/ such that 

i = i + 27rx2 S . (12) 

where s is an arbitrary integer. There, we showed that the physical observables 
have exactly the same matrix elements in both theories and, consequently, that 
the theories are physically equivalent. In particular, a theory of interacting 
fermions is always equivalent to a family of theories of interacting fermions 
coupled to a Chern-Simons gauge field with coupling constant 9 such that 
a = 2tt x 2s. This result is the starting point of our analysis of the FQHE. 

From Eq it is apparent that the physics of these theories is a periodic 
function of the statistical angle, <5 — > S + 2tt x integer. Even though these the- 
ories are exactly periodic, periodicity is broken in approximate, semiclassical 
approximations, such as in the Average Field Approximation that we will use 
here since one particular period will have to be chosen. This is an important 
issue since, at the perturbative level, the fluctuations around this mean-field 
are not able to restore periodicity as an exact property. The restoration of 
periodicity is a non-perturbative effect. In particular, operators which create 
excitations that change the amount of flux per particle in even multiples of 
the flux quantum are soliton operators which restore the periodicity broken 
by the mean-field. However, we can take advantage of the periodicity of the 
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CS description to choose the period in which the mean-field-theory is simplest. 
This is the approach we take in this work to attack the FQHE. As it will be 
apparent in the next Section, at the mean-field level, this approach reproduces 
Jain's construction of the FQHE states. Recently, Kivelson, Lee and ZhangEj 
have forcefully argued that periodicity plays a central role in the determination 
of the global phase diagram of the 2DEG in a strong magnetic field. 

In second quantized language, what we have proven is the physical equiv- 
alence of theories whose dynamics are governed by the actions ( in units in 
which h = 1) Sg and Sgi which are defined by 

S e = Jd 3 z ^%z)[iD +M(z)-^\D4>(z)\ 2 + ^ x a*F x \ 

- ljd 3 zj d 3 z' (|^(z)| 2 - p)V(\z- ?\)(W)\ 2 ~ P)- (13) 

where p is the average particle density, provided that 9 and 9' satisfy Eq. |lj. In 
Eq. ^ V'( 2; ) is a second quantized Fermi field, /i is the chemical potential and 
.Dp is the covariant derivative which couples the fermions to both the external 
electromagnetic field and to the statistical gauge field o M 

D /Ji = d fl +i^A fi +ia fi . (14) 

In particular, a theory of interacting fermions ( which has s — 0) is equivalent 
to a family of theories of fermions with \ = 2irx even integer. 



2.1 The Semiclassical Limit and The Laughlin Ground State 

In the remaining of this Section we will show that the semiclassical limit of the 
theory described by the action Se of Eq.|l^, with \ — 2ir(m— 1), yields the same 
physics as the Laughlin state. In order to prove this statement we will develop 
a semiclassical approach to this problem. As a by-product, our formalism 
provides for a systematic procedure to compute corrections to the Laughlin 
approximation. This is, to the best of our knowledge, the first formalism for 
which the Laughlin ansatz arises as the first of a series of approximations. 
Consider the quantum partition function for this problem ( at T = 0) 

Z= vrv^Va^MiSe). (15) 

We will treat this path-integral in the semiclassical approximation. In order to 
do that, we will first integrate-out the fermions and treat the resulting bosonic 
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theory within the Saddle Point Expansion (SPE) characteristicxifjSemiclassical 
approaches to quantum mechanics and quantum field theory liilrLa. In terms 
of Feynman diagrams, the semiclassical expansion is an expansion in powers 
of the number of fermion loops (bubbles). Thus, the leading semiclassical 
approximation (gaussian fluctuations) is equivalent to the conventional RPA 
of many-body physics. It is easy to see that the formal expansion parameter 
that organizes the SPE is the filling fraction v. This is not surprising since 
RPA is a high density approximation E3. However, for the FQHE, the filling 
fraction v is never large but typically of order one. We will show below that, 
in spite of this, the SPE will still be very useful since, for special values of the 
filling fraction, it will enable us to construct a set of ground states for which the 
spectrum is fully gapped. Thus, for those states the SPE yields a perturbative 
expansion free of infrared divergencies. However we will also encounter a set 
of fractions for which there is no gap and the SPE (or equivalently, RPA) is 
infrared divergent at every order. Interestingly enough, the set of fractions for 
which the SPE predicts a fully gapped state coincides with the main fractions 
of the FQHE. These are the states in Jain's hierarchy. 

In the absence of electron-electron interactions, the fermions can be in- 
tegrated out immediately since the action becomes quadratic in Fermi fields. 
In the presence of interactions this is no longer possible since the interaction 
term in the action Eq. [l3] spoils this feature. However, the physical states 
used in deriving the path integral must obey the local charge-flux constraint 
of the Chern-Simons theory. Hence, it is legitimate to replace the charge den- 
sity jo(x) by 9B(x) in the pair-interaction term of the action, at all points of 



space-time x and to write the pair-interaction term of Eq. 13 in the form 



S in t = \f d3z J d 3 z'(8B(z) - p)V(z - z')(6B(z') - p)) (16) 

where V(z — z') represents the instantaneous pair interaction i. e., 

V(z-z') = V(\z-?\)5(t-t'). (17) 

We will assume that the pair potential has either the Coulomb form, i. e., 
^(H) = V> or that in momentum space it satisfies that V(Q)Q 2 vanishes at 
zero momentum. This includes the case of ultralocal potentials, (i. e., with a 
range smaller or of the same order as the cyclotron length £), in which case we 
can set V — 0, or short range potentials with a range longer that £ such as a 
Yukawa interaction. 

Note that since we are dealing with a gauge theory, a gauge has to be 
specified in order to make the functional integral well defined. We will assume 
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from now on that a gauge fixing condition has been imposed but, for the 
moment, we will not make any specific choice of gauge. 

The partition function Z can be written in the form of a functional integral 
involving the Fermi fields ip, and the statistical gauge fields a M . The Fermi 
fields can be integrated out without any difficulty yielding, as usual, a fermion 
determinant. The resulting partition function can thus be written in terms of 
an effective action S c fi given by 



Simons action Eq. |j. The field represents a small fluctuating electromag- 
netic field, with vanishing average everywhere, which will be used to probe the 
system. The electromagnetic currents will be calculated as first derivatives of 
Z with respect to A^. The full electromagnetic response will be obtained in 
this way. Notice that we have used the invariance of the measure T>a^ with 
respect to shifts, to move A^ out of the covariant derivatives and into the 
Chern-Simons term Scs- 

We are now ready to proceed with the semiclassical approximation. The 
Saddle Point Equations, or classical equations of motion are 



By varying S c g with respect to we find 



where {jfi(z))p represents the expectation value of the charge and current of 
the equivalent fermion problem. 

These equations have many possible solutions which include uniform (liq- 
uid) states, Wigner crystals, and non-uniform states with vortex-like con- 
figurations. We will only consider solutions with uniform particle density 
{{jo{z)) — p), i.e., the liquid phase solution, and no currents in the ground 
state. 




(M z ))f 
{jk(z)) F 
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If the external electromagnetic fluctuation is assumed to have zero average, 
the only possible solutions of this type are 

w - -| 

(£) = (21) 

This is the average field approximation (AFA) which can be regarded as a mean 
field approximation. Eq. |2l] shows that, for a translationally invariant ground 
state, the effect of the statistical gauge fields, at the level of the saddle-point- 
approximation, is to change the effective flux experienced by the fermions. 
The total effective field is thus reduced from the value of the external field B 
down to £? c ff = B + (B) = B — |. We want to find the ground state and the 
spectrum of low energy excitations of a system of N (interacting) electrons 
in the presence of an external magnetic field of strength B. We will further 
assume that the linear size L of the sample is such that a total of quanta 
of the magnetic flux are piercing the surface. In general, the filling fraction 
v = jf- is not an integer. Notice that, as B and p are varied, the effective field 
-B e ff may either be parallel or anti-parallel to the external field B. Thus, we 
will not assume a particular sign for B c g although we will set B > 0. 

The uniform effective magnetic field B c g, which solves Eq. |2l], define a 
new set of effective Landau levels. Each level has a degeneracy equal to the 
total number of effective flux quanta N e g and the separation between levels is 
the effective cyclotron frequency — e ^ c " . Similarly, there is an effective 
cyclotron radius £ eS . It is easy to see that the uniform saddle-point state 
which satisfies Eq. |2l] has a gap only if the effective field B e g experienced by 
the N fermions is such that the fermions fill exactly an integer number p of the 
effective Landau levels. This is precisely the point of view advocated by Jain: 
the FQHE is an IQHE of a system of electrons dressed by an even number 
of flux quanta. However, this condition cannot be met for arbitrary values of 
the filling fraction v at fixed field ( or at fixed density). Let N^ s denote the 
effective number of flux quanta piercing the surface after screening. It is given 

b y 

±2nNf = 2ttA -|l 2 . (22) 

where the ± sign stands for the case of an effective field parallel or antiparallel 
to B. Thus, the effective cyclotron frequency is reduced from its free 
electron value of down to wf = uj c (1 — j^g)- The effective cyclotron 
radius is given by £ off = (£/ '«/! — j^g) which is larger than the non- interacting 
value. Therefore, even though the bare Landau levels may be separated by a 
sizable Landau gap fw c , the effective Landau levels have the smaller gap fru^ , 



10 



Substituting the value of 9 in Eq. we obtain 

± 2nNf = 2vrA^ - 2n2sN. (23) 

where 2s is an even integer. The spectrum supported by this state has an 
energy gap if the N fermions fill exactly p of the Landau levels created by the 
effective field B c s . In other words, the effective filling fraction is f e ff = tt^tt — 

p. Using Eq. 03], we find that the filling fraction v and the external magnetic 
field B must satisfy 

N N 

± — = 2sN, (24) 

p v 

or, equivalentlythe allowed filling fractions are 

^ = 2^L1 (25) 

The allowed filling fractions v± can thus be arranged into families or hi- 
erarchies. As p — > oo, the allowed fractions approach the limiting values 
lmip^oo v±(j3, s) = i-, with v + {p, s) < while v-(p, s) > j^. For the special 
case s = 1, the limiting value is |. Since v+{p, 1) < h and 1) > i, the 

sequence v+(p, 1) is said to represent electron-like FQH states whereas the mir- 
ror sequence V-(p, 1) is said to represent hole-like FQH states. The effective 
Landau gap for these solutions is 

fi »" = ^h < 26 » 

which is small if either p or s are large. Thus, the energy to excite a fermion 
is Ulo^ s and it is considerably smaller than the bare free-particle value hu) c . 

The states are thus parametrized by two integers p ( the number of filled 
effective Landau levels of the effective field) and 2s ( the number of flux quanta 
carried by each fermion), and by a sign. It is trivial to see that the Laughlin 
sequence is an obvious solution of Eq. [2^ for p = 1 and 2 s = m— 1 . The effective 
fermions thus fill up exactly one Landau level and has to be chosen to be 
\ = 2ix(m — 1). This is Jain's result. At this mean-filled level the wave function 
is the Slater determinant for one filled Landau level xi of Eq. |^. We will show 
below that the additional factor W^Azi — zf) m ~ x is due to fluctuations. For 
general p and s, these are the allowed states in Jain's hierarchy uu. Please 
notice that, unless s = 0, all states in these hierarchies have v±(p, s) < 1. For 
s = (i. e., no flux attachement) v = p and we reproduce the integer QH 
states. Thus, only primary FQH states are generated. States characterized 
by other fractions can be constructed by means of the conventional scheme of 
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condensing quasiparticles or quasiholes on top of these scqucncics. We will see 
below that, as expected, the quasiparticles and quasiholes of the FQH states in 
these hierarchies are anyons. However, it is unclear how many of these higher 
level hierarchical states are energetically stable (or accessible). 

The saddle-point approximation yields a very simple spectrum which con- 
sists just in the single and many-particle excitations of fermions in the effective 
field B e g. The single-particle gap is equal to Huj^ = hu c /{2sp ± 1). Thus, 
at fixed B, the single particle energy gaps become smaller as the index p in- 
creases and vanish as p — > oo. In this limit the hierarchical states converge to 
the fractions v — For these fractions the system is compressible, i. e., there 
is no energy gap in its spectrum, and the perturbation expansion breaks down. 
The breakdown is signalled by infrared divergent contributions at every order 
of the expansion at low tenaaeratures. The case of the ^ states was analyzed 
by Halperin, Lee and Readta. Elsewhere in this volume, Steve Simon presents 
a detailed discussion of this important case. 

At this level of the SPE, the Coulomb energy does not enter explicitly in the 
effective single particle gaps. However, their existence is a physical consequence 
of the existence of Coulomb interactions. Moreover, in the strong field limit 
B — ► oo, both the bare and the effective Landau gaps diverge. This is clearly 
incorrect since, as B increases the states of the system should be made almost 
completely of states only in the lowest Landau level. In the strong field limit, 
the Coulomb energy alone should determine the numerical value of the gaps. 
The origin of this difficulty is that the SPE (or the RPA) is formally correct 
in the small field limit and the states exhibiting the FQHE are all in the large 
field limit. Thus, we should expect that the higher order correction in the SPE 
should not only correct the numerical value of the gaps but also lead to a finite 
limit in the strong field regime. These corrections have not been calculated 
and their properties are not yet understood. The same cons ider ations apply 
to the spectrum of collective excitations discussed in Section 3J. However, it 
is important to stress that the SPE (with gaussian corrections) yields a set 
of states with the correct quantum numbers even though the energies are not 
correct. In particular we will show below that, once gaussian fluctuation are 
accounted for, these single particle states have the correct charge and statistics 
(although the wrong energy). 

The saddle-point equations, Eq. |2^, have a host of non-uniform solutions 
which have finite energy. They can be viewed as soliton or vortex solutions. 
Our construction is very close in spirit to the picture of the quasihole presented 
by LaughlinllZl. We will only discuss here the qualitative features of these 
solutions. For the sake of simplicity, we will only study the case of fermions 
in the Laughlin — sequence. Recall that this sequence is represented by the 
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solution with p — 1 and m = 2s + 1. The uniform state was constructed 
by filling up the lowest effective Landau level. Let us consider the state which 
results from removing a fermion from the single-particle state, centered around 
the origin z = and with lowest angular momentum, and placing it on the 
first unoccupied angular momentum state. Physically, this new state lies on 
the outer edge of the system. For an uniform effective field, this state does not 
exist. But, if the effective field is increased at the origin by an amount equal 
to one flux quantum, the angular momentum of all its eigenstates is raised by 
one whole unit. The radius Rn of the droplet with N particles swells to a new 
value Rn + SR large enough to include a new cyclotron orbit. Qualitatively, a 
quasihole localized at zq has the mean-field wave function = ^!h(zo, {zj}) 
given by 

N 

* h (z , { Zj }) = Y[( Zi - z ) Xi({zj}). (27) 

i=l 

where Xi({ z j}) is the wave function for TV fermions occupying the lowest Lan- 
dau level. Notice that this wave function differs from the Laughlin state for 
the quasihole by the prefactor Ili^' ( z » — z j) m ~ 1 - Indeed, this prefactor is also 
missing in the mean-field wave function for the ground state. In our picture, 
both prefactors arise from fluctuations which attach fluxes to the particles. 
The quasihole wave function of Eq. is an approximation valid in the limit 
\zi — zq\ 3> i cS - From Eq. it is easy to see that the mean- field excita- 
tion energy of the quasihole Eh is given approximately by the Coulomb energy 
V{l cS ). 



2.2 The role of the Gaussian fluctuations 

We consider now the role of the gaussian fluctuations around the uniform clas- 



sical solutions discussed in Section 2.1. We begin by considering the effective 
action of Eq. Let a M (x) denote the fluctuations of the statistical vector 
potential a^, i. e., we set a M — > (a M ) +d^. The effective action of Eq. fl8| can be 
expanded in a series in powers of the fluctuations. We will be interested only 
in keeping just up to quadratic terms in the fluctuations. In the language of 
Feynman diagrams, we are summing up all the one-loop bubble contributions. 
Thus, the effective action at the quadratic level involves the linear response 
kernels (evaluated in the RPA) for a system of fermions in an external static 
uniform magnetic field B c g with an integer number of filled effective Landau 
levels. As usual, the linear terms are cancelled if the saddle-point equations 
are satisfied. 
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At the quadratic (gaussian) level the effective action has the form 



^ (2) = \j ^AM*) H^foy) a„(y)+5cs(a M -i„)+C 3 (^-4)- ( 2 ») 
where 

Slnt (2,. -Ay) = ~ J d 3 xd 3 y [b(x) ~ B(x)] V(x - y) [%) - B{y)\ (29) 

Notice that the only approximation used in Eq. is in the first term which 
follows from the expansion of the fermion determinant in powers of a^. The 
remaining terms are exact. 

Eq. ^ holds provided that the non-quadratic dependence in the fluctu- 
ating part of the statistical vector potentials a M is small. Recall that these 
non-quadratic terms result from expanding the (logarithm) of the fermion de- 
terminant in powers of the fluctuations around the average field approxima- 
tion. The kernels that enter in the expressions for these terms are (connected) 
current correlation functions (or response functions) of the mean field the- 
ory. Thus, the tensor H ®} is the polarization tensor of the equivalent fermion 
problem at the mean field level, and it is obtained by expanding the fermion 
determinant up to quadratic order in the statistical gauge field 

Z [A„] = Det [iD + (i + A + ^D 2 } 

= Z[0] exp{^ J d D x J d D y A^x)Il^(x,y)A u (y) + . . .} (30) 

The tensor Ilffl should not be confused with the true electromagnetic polar- 
ization tensor which measures the response of the whole system to a weak 
electromagnetic field. We will come back to this issue in the next Section. 

It is straightforward to derive an expression for n<3 in terms of the one- 
particle Green functions G(x, y) 

G(x,y)=<x\— ] - \y>. (31) 

iD + fi + ^jD 2 [< A >} 

The components of the polarization tensor Ilffi(x,y) are (ft, = 1) 

Hoo (»>?/)= iG(x,y)G(y,x) 

<W) = ^{G{x,y)D)G{y,x)-G{y,x)DfG{x,y)} 
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IlfoW) = +^{~G(x,y)DfG(y,x)+G(y,x)D*G(x,y)} 
rf$(*,v)= ^S(x-y)6 jk G(x,y) + 

i 



~ ^( D ! G (^y))( D t G (y^)) - j^{DfG(y,x))(DfG(x,y)) + 

+ ^2 G (V' x))(D]D^G(x, y)) + D V k G(y, x))G(x, y) ■ 

(32) 

It can be shown that the effective action of Eq. |3^ is gauge-invariant and that, 
in consequence, Tlffi is transverse, i. e., 

d*ll$(x,y)=0. (33) 

In reference we give detailed expressions for the various components of IT^ J . 
In particular we show that the action is gauge-invariant and that IljS is trans- 
verse, albeit only in a weak sense i. e., as a distributkan. Eq. ( |32| ) agrees with 
the results of Randjbar-Daemi, Salam and StrathdeeE3 but disagrees with the 
calculation of II^ J by Chen et. al. El. 

From the expressions for the components of nJS in referenced, we see that 
this tensor is analytic in -g-^ and that it has simple poles at ui — ku c (with k 
an integer), where Co c = uj c /(2sp+ 1) is the cyclotron frequency associated with 
the effective magnetic field -B e ff- As a result, II^J has a gradient expansion 
in powers of the inverse of the effective magnetic field -g^— , or equivalently, in 

powers of the inverse of the external magnetic field -g. In fact, the dimension- 

o 2 

less parameter of this expansion is (we are working in a system of units 
such that h = c = e = 1). It also turns out that, within this approximation, 
the limits of B — > oo and M — > are not equivalent (see the explicit form of 
njS given in referencei). 

The non-quadratic terms in in the effective action are of the form 

S cS = S^ia-^A^ + ^ : Jd 3 x d 3 y d 3 z S" ~a v a x I$ x (x,y,z) + ... (34) 

where the kernel Tl^ x (x, y, z) represents a three-point current correlation func- 
tion in the mean field theory. Thus, in the language of Feynman diagrams, 
while n)2 (x, y) can be viewed as a fermion bubble with two amputated ex- 
ternal collective mode lines, H^ x (x,y, z) again has one fermion loop tied to 
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three amputated external collective mode lines a M . Each one of these non- 
quadratic kernels have the same gauge invariance (i. e., transversality) and 
analytic properties as the gaussian (or RPA) kernel. In particular, this means 
that, in momentum and frequency space, these kernels must be a linear combi- 
nation of tensors (of the appropriate rank) which have the correct transversality 
properties, times a set of functions which are analytic in Q 2 and have poles 
at frequencies equal to an integer multiple of the effective cyclotron frequency. 

Therefore, the non-quadratic terms necessarily have powers of -g— (for each 
one of the external momenta and frequency entering the fermion loop) which 
are higher than the ones found at the quadratic level. Since the mean field 
theory has an integer number of filled Landau levels, the energy denominators 

of the kernels do not change this counting in powers of In conclusion, the 
expansion of the fermion determinant, and hence of the effective action, is actu- 
ally an expansion in powers of -§—, or equivalently, in powers of However, 

an expansion in powers of is also a gradient expansion. Thus, the gradient 
expansion and the semiclassical expansion mix and are not independent from 
each other. 

The semiclassical expansion is obtained according to the following rules. 
The propagator for the fluctuations, which represent collective modes, is the 
inverse of the kernel of the gaussian action. Since the pair potential enters 
only through the propagator for the fluctuations, the perturbation theory is 
not an expansion in the powers of the pair interaction. From this point of 
view, this expansion is very different from conventional expansions around 
the Hartree and Hartree-Fock approximations. The vertices of the expansion 
are the kernels for the non-quadratic terms. This expansion lacks a natural 
small parameter (i. e., a coupling constant) and it should be regarded, like 
all semiclassical expansions, as an expansion in the number of fermion loops 
(i. e., RPA plus corrections). One should keep in mind, however, our previous 

discussion on its exactness in powers of ^g-. In what follows we will make 
extensive use of the formal properties of this expansion. 

We remark here that, unlike other mean field approaches to conventional 
many body systems (such as Hartree-Fock), the gaussian corrections of this 
theory must alter the qualitative properties of the state described by the AFA. 
The reason is that the AFA violates explicitly space-time symmetries, such as 
Galilean invariance (more generally, magnetic invariance) which, for transla- 
tionally invariant systems, must remain unbroken and unchanged. Thus the 
center of mass of the system must execute a cyclotron- like motion at, exactly, 
the cyclotron frequency of non interacting electrons in the full external mag- 
netic field, as demanded by Kohn's theoremEj. A naive application of the AFA 
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would suggest that the cyclotron frequency is renormalizcd downwards since 
the effective field seen by the composite fermions is smaller than the external 
field B. Hence, the magnetic algebra may appear to have changed. We will 
see below that the gaussian fluctuations yield the correct cyclotron frequency 
and, thus, restore the correct magnetic algebra. 

3 Response functions and spectrum of collective excitations 

The physical systems which exhibit the FQHE present a very rich response to 
external electromagnetic perturbations. While some of the observed phenom- 
ena, such as cyclotron resonance, can be understood in terms of simple global 
motions of the center of mass under the combined influence of electric and 
magnetic fields, the spectrum of collective excitations is certainly determined 
by the interactions. Given the unusual features of the Laughlin states and its 
generalizations, it is expected that the same features should largely determine 
the behavior of the collective modes too. We will now use the fermion field 
theory described in Section ||, to study the collective excitations of the fully 
polarized FQHE states in the sequence 1/v = ±l/p + 2s (where p and s are 
two positive integers). While the average field approximation is certainly very 
appealing, it has the serious problem that it breaks a number of space-time 
symmetries quite explicitly. In particular, it breaks both Galilean and Mag- 
netic invariance. It turns out that the leading quantum fluctuations around 
this state, restore these symmetries, in the uniform Q — > limit, already at 
the gaussian level. Indeed, we find that the quadratic or gaussian level of the 
semiclassical expansion gives the correct value of the Hall conductance of the 
system, and that the leading order of the density correlation function satu- 
rates the /-sum rule. This is an essential result to show that the absolute value 
squared of the wave function of all the (incompressible) liquid states has the 
Laughlin form at very long distances, in the thermodynamic limit (see Section 



3.1 Electromagnetic response junctions for the FQHE 

The effective action for the fluctuations of the Chern-Simons gauge field (a M ), 
at the Gaussian level of our approximation, is given by 



D- 



2 



1 



/ 



d 3 x d 3 y a»(x) IL$(x,y) ~a v {y) + -S cs (^ - A„) 




d 3 x d 3 y B(x) - B(x) V(x - y) B{y) - B(y) 



(35) 
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Since is quadratic in a M , we can integrate out this field and obtain the 
effective action for the electromagnetic fluctuations A^, S^g(A /1 ) . We will use 
this effective action to calculate the full electromagnetic response functions at 
the gaussian level. Since this calculation is based on a one loop effective action 
for the fermions (i. e., a sum of fermion bubble diagrams), this approximation 
amounts to a random phase correction to the average field approximation. 

In order to integrate out the statistical gauge field a M we must fix the gauge. 
The electromagnetic effective action, being gauge invariant, is independent of 
the choice of gauge for the statistical gauge fields in the path integral. We fix 
the gauge d^a^ — a (where a is an arbitrary real number) using the standard 
Faddeev-Popov procedure []. The result is explicitly gauge invariant and all 
dependence on the parameter a cancels out. At the one-loop level (governed by 
the effective action of Eq. |35| ) we need to know the inverse of the polarization 
tensor of the equivalent fermion problem, IT^. In referenceB, we showed that 
IljS can be written in terms of three gauge invariant tensors, an E 2 term, a B 2 
term, and a Chern-Simons term. These three tensors plus BV.E and a gauge 
fixing term (such as ^(c^a^) 2 which corresponds to the Landau-Lorentz gauge 
if a — > 0), close an algebra that can be used to invert the polarization tensor 
and to calculate explicitly the electromagnetic response functions. 

After integrating out the statistical gauge field in Eq. the effective 
action for the electromagnetic fluctuations turns out to be 

«5 e T(4) = \J d 3 xd 3 yA^x)K^(x,y)A u (y) (36) 

Here K^ v is the electromagnetic polarization tensor. It measures the linear 
response of the system to a weak electromagnetic perturbation. Its components 
can be written in momentum space as follows 

^oo = Q 2 K (u,Q) 

K 0j = LoQjK (u),Q) +ie jk QkKi(Lu,Q) 

Kj = loQ 3 K (u, Q) - iejkQkK^oj, Q) 

K i3 = u 2 6ijK (u),Q) - ieijU)Ki(u,Q) + (Q 2 6ij — QiQj)K 2 (uJ,Q) 

(37) 



a To be more precise, whenever we have had to fix the gauge we have adopted Feynman's 
method of averaging over gauges, with a being the width of the distributions. Thpiis a 
standard procedure which is reviewed in many textbooks (see, for instance, reference cj) • 
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where the functions Ki(ui, Q) (I = 0, 1, 2) are given by 



K (u, Q) 



n 



D(u,Q) 

2 (fl+ni) , rtrr,KA2 



Ki(u,q) = e + e z \; ^ +e 6 v(Q)Q 



D(u>,Q) D(u,Q) 

n 2 v(Q) (lu 2 n 2 - n! 2 ) q 2 n n 2 



and 



(38) 

D(w, Q) - n V - (ni + ef + n (n 2 - e 2 v(Q))Q 2 (39) 



The coefficients II; (/ = 0, 1, 2) are functions of cj and Q, and are given explicitly 
in referenceS. V(Q) is the Fourier transform of the interparticle pair potential. 
As we mentioned before, we needed to include a gauge fixing term to be able 
to compute the functional integral in Eq. |35|. But at the end of the calculation 
all the terms which contain the gauge fixing coefficient (a) cancel each other, 
and the final result for the response functions is, as it must be, gauge invariant. 
The other tensor that we have introduced to make the calculations, BV.E, is 
not present in the final answer either. 

We want to stress here that the thermodynamic limit is crucial for the 
accuracy of our results. Notice first that in the electromagnetic effective action 
of Eq. |3^ we are neglecting higher order response functions, i. e., correlation 
functions of three or more currents or densities. We have shown in Section 



2.2 that these higher order correlation functions have higher order powers of 



^fc- than the quadratic term. Strictly speaking, these terms are not neglectible 
for a finite system because, in this case, there is a minimum value that the 
momentum can take, determined by the linear size of the system L, i. e., 
\Q\ > But in the thermodynamic limit, L — > oo and the momentum can go 
to zero. In other words, only for an infinite system one is allowed to keep only 
the quadratic term in the electromagnetic action, Eq. |36|, and to neglect the 
higher order correlation functions. 

The electromagnetic response functions determined by K^ v have the fol- 
lowing properties: 

i) The polarization tensor at mean field level, Hj^J , has poles at every 
value of the effective cyclotron frequency (ui c = ^ff-)- This corresponds to the 
physical picture, at mean field level, of an IQHE of the bound states in the 
presence of a partially screened external magnetic field, B c g . Once we take into 
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account the Gaussian fluctuations, it is easy to prove that all this poles that are 
present in the numerator and the denominator of the K^y components through 
the Ilj's, cancel out, and the poles of the response functions are determined 
only by the zeroes of their denominator, D(ui, Q). In other words, the collective 
excitations of this system will be determined only by the zeroes of D(ui, Q). 

ii) The leading order term in Q 2 of the Kqq component of the polarization 
tensor saturates the /-sum rule. 

iii) The gaussian fluctuations of the statistical gauge field are responsible 
for the FQHE. In particular, the gaussian corrections yield the exact value for 
the Hall conductance. 

In the remaining of this section we will discuss these properties and their 
experimentally accessible consequences in detail. 



3.2 The Spectrum of Collective Excitations 

For simplicity, we give here the zeroes of D(uj, Q) for the case of p — 1 filled 
effective Landau levels. The (more tedious) case oLgeneral p can be studied 
by straightforward application of the same methodsEj. 

In the case p = 1 the filling fraction is v = where m = 1 + 2s, i. e., the 
Laughlin sequence. We find that there is a family of collective modes whose 
zero-momentum gap is kui c , where k is an integer number different from 1 and 
m, and whose dispersion curve u>k(Q) is 



\k-l -,2 



2fc(m-l)(fc-l) 



'2B eS ' c (k - 1)! (k - m) 
The residue in Kqq corresponding to this pole is 



Res(K 00 ,cj k {Q)) = -Q 2 uj c ^- (-§-) 



\fc-l 



2fc(m-l)(fc-l) 



(40) 



(41) 



'2B eS ' (fc - l)!(fc - m){k 2 - n„ , 

The cases k = 1, m have to be treated separately. In general, we find that 
there is no mode with a zero momentum gap at Q c . Instead, at Q = 0, there 
is a doubly degenerate mode with a gap at u> c . This degenerate cyclotron 
mode can be viewed as the mixing of the modes with k = 1 and with k = m. 
Thus, the mode with k = 1 has been "pushed up" to the cyclotron frequency 
(at Q = 0). For Q ^ 0, the degeneracy is lifted and these two modes have 
different dispersion curves. 

For the special case of v — -|, i. e., m = 3, this effect is particularly 
important. The dispersion relations for the cyclotron modes are given by 



w±(Q) 



^ 2£> c ff ^ 



) T oy= 



(42) 
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where 



a± 



2MV(0) 
2^ 



± 



2MV(0) 2 



2tt 



The residues corresponding to these poles are 



Res(K 00 ,u±(Q)) 



2tt 



(1 



Y + 288 



288 s 



(43) 



(44) 



For v = to > 5, the collective modes whose zero- momentum gap is the 
cyclotron frequency, u> c , are 



LU+(Q) 

with residue 



-2 ( TO 

x 4 V 



(to - 2) 



2A/V(0)\ 
2^ / 



ites(if oo, 



-Q Wc - 



(45) 
(46) 



The other cyclotron mode has the dispersion 



w_(Q) 



(-^-) 



m— 2 -2 



4m 2 (m-l) 2 /. (m-1) 2MV(0)y 



(to - 1)! 



(m - 2) 



2tt / 



(47) 



with residue 
Res{K 00 , w_(Q)) 

rj2 



— ( 



2tt v 2B cff 



3 4to 2 (to-1) 2 f 4 (m- 1) 



(to - 1)! 



(m-2) 



2^ J 

(48) 



The above results are valid only if the pair potential V(Q), has a gradient 
expansion in powers of Q, i. e., for short range interactions. V(0) stands for 
the leading order term in that expansion. 

If the pair potential has the Coulomb form, i. e., V(Q) = in two 

spacial dimensions, both, the dispersion relations with zero-momentum gap at 
the cyclotron frequency and their residues get modified. The expressions valid 
in this case are, for any allowed value of m 



w+(Q) = 



\Q\ 

2B cS 



Uj\ 2Mq 2 



with the same residue given by Eq. Em and 



w_(Q) 



\Q 



2m-3 



Am 2 (m - l) 2 
(2B eS ) m - 2 ~ c 2Mq 2 (m- 1)! 
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(49) 



(50) 



with residue 



igi 2 (" 1 - 2 ) 4m 2 (m-l) 2 
Res{K 00 ,W-{Q)) = -Q 2 . c - x ^ 2 Mg 2 ) 2 (m - 1)! ^ ! 



In this section we have found the spectrum of collective excitations for 
some values of the fillineJjraction. Our results are a generalization of the work 
of Kallin and HalperincJ who studied the spectrum of collective modes for 
the integer quantum Hall effect within the RPA. We find a family of collective 
modes with dispersion relations whose zero- momentum gap is ku) Cl where k 
is an integer number different from 1 and m. When k = m, i. e., the zero- 
momentum gap is the cyclotron frequency, there is a splitting in the dispersion 
relation for finite wavevector. One of these solutions, has negative slope 
for small values of Q. Therefore, there must be a roton minimum at some finite 
value of the wavevector. Since our results are accurate only for small Q, our 
dispersion curves do not apply close to the roton minimum. Nevertheless, this 
mode is expected to become damped due to non-quadratic interactions among 
the collective modes. On the other hand, the collective mode with lowest 
energy, which has k = 2, is stable (at least for reasonably small wavevectors) 
and, at small wavectors, it disperses downwards in energy. This behavior 
suggests that there should be a magnetoroton minimum for this mode. This 
result is consistent with the work of Girvin et al. E3, where the authors used the 
single mode approximation to obtain the lowest collective mode in the lowest 
Landau level (the intra mode) for the states in the Laughlin sequence. The 
generalization of these .results to other states of Jain's hierarchy (p ^ 1) has 
been given in referenceE 2 ! Our results coincide, in the small momentum limit, 
with the ones obtained by Simon et al. E3 more recently. The splitting of the 
cyclotron mode for v = 1/3 is a little puzzling. It only happens for v — 1/3 and 
for short range interactions. In all other cases, only the residue for one of the 
two cyclotron modes is proportional to Q 2 . Standard lore has it that Kohn's 
theorem demands that there should be one and only one mode converging to 
the cyclotron frequency as Q 2 =± with residue proportional to Q 2 . Zhang has 
emphasized this point recentlylHj. It is generally assumed that Kohn's theorem 
is valid even at non-zero wavevectors and that it requires the existence of only 
one mode with residue proportional to Q 2 converging to lo c . However, at 
non-zero wave vectors, these arguments make the unstated assumption of the 
analyticity of the current operators on the wavevectors. While this may well be 
correct, it is an additional assumption and it does deserve closer scrutiny. The 
results from our theory do indeed predict the existence of only one mode at 
ui c with residue proportional to Q 2 , which is the statement of Kohn's theorem. 
And, also, for all filling fractions and for all pair potentials (except v = 1/3 
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and short range interactions) we do find only one mode with residue Q 2 even at 
non-zero wavevectors. The case v = 1/3 and short range interactions appears 
to be exceptional in that we find two modes which coalesce at the cyclotron 
frequency as Q 2 — > 0. But both of these modes have residue proportional to 
Q 2 , with different amplitude, and together they satisfy the sum ruleEJ. While 
it is possible that the non-gaussian corrections may change this result since, in 
a sense, these are subleading pieces in Q 2 , these non-gaussian corrections are 
expected to be very small at small wavevectors. 

We close this section with a few comments on the validity of this spectrum 
of collective modes beyond the semiclassical approximation. Primarily we have 
to consider the physics at moderately large wavevectors and the (expected) 
effects of non-gaussian corrections. At the gaussian (RPA) level we found a 
family of collective modes which, for sufficiently small momentum, arc infinitely 
long lived (i. e., the response functions have delta- function sharp poles at their 
location). These modes represent charge-neutral bound states. It is in principle 
possible that, for Q sufficiently large, these modes should become damped. The 
threshold should occur when the energy of the collective mode becomes equal 
to the energy necessary to create the lowest available two-particle state: a 
quasiparticle-quasihole pair. In the AFA, the energy of a pair is equal to ui c . 
Gaussian fluctuations are expected to renormalize this energy and to give it 
a momentum dependence. This is in principle calculable with the methods 
described here El Non-gaussian corrections to the RPA are also expected to 
give a finite width to (presumably) all the collective modes but the lowest one. 
This is so because the corrections to the semiclassical approximation are due to 
effective vertices (due to virtual quasiparticle-quasihole pairs) which couple the 
various collective modes and, thus, induce the higher energy modes to decay 
down into the lower modes. However, by gauge invariance, these vertices have 
a momentum dependence and should vanish as Q — > 0. Thus, the width of 
the higher energy modes goes to zero as Q — > and these modes only become 
sharp at Q = 0. But at Q = the only accessible mode is the cyclotron mode 
(the other modes have a vanishingly small spectral weight). These arguments 
strongly suggest that the only truly sharp mode, at Q =f~-P, is the cyclotron 
mode, which is required to be stable by Kohn's thcorcmE3. Since the modes 
with zero momentum gap at kuj c , k > 3, are not the collective modes with 
lowest energy, it is possible that at finite wavevectors they may also decay into 
the collective mode with lowest energy ( the mode with k = 2, which has a 
gap at u c ). 
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3.3 Optical properties and experiments 

In this section we discuss the experimental consequences of the results that we 
have just derived. The density correlation function can be probed by optical 
absorption and by Raman-scattering experiments. 

In the first case, the optical absorption is proportional to the imaginary 
part of the density correlation function. We predict that there will be ab- 
sorbtion peaks at a discrete set of frequencies which, for Q — > converge to 
lu = fc(D c , where k is an integer number greater than two. Since the spectral 
weight of these modes vanishes as Q — > 0, the associated absorption peaks 
are, for a strictly translationally invariant system, only observable at non-zero 
momentum. 

In the case of the Raman scattering, the geometry must be such that there 
is a component of the incident light wavevector in the plane of the sample. 
The Raman spectrum, /(w) ( _is also proportional to the imaginary part of the 
density correlation function E3. 

We have seen that in the limit \Q\ <C where I is the magnetic length, 
most of the weight of K 00 (lu,Q) is in one of the cyclotron modes. The pole 
in Kqq(u!, Q) for the lowest excitation frequency, lo^ with k = 2, has a residue 
which is proportional to \Q\ A , i. e., it is smaller by a factor of \Q\ 2 than the 
residue at the highest weighted mode at the cyclotron frequency. 

We have also found that there is a splitting in the cyclotron modes. If the 
pair potential has a gradient expansion in \Q\, i. e., short range interaction, the 
pole at id- (Eq. |47]), has a residue that is smaller by a factor of |Q| 2( - m_3 ^ than 
the residue of lo + (Eq. ^5|) . The relative Raman intensity, i is proportional 
to (25fi)(™- 3 ) which is a big number within our approximation. If the filling 

fraction is v — i, both modes have the same Q 2 dependence in their spectral 
weight, but the relative intensity is w 2.5 provided that V(0) = 0. Except for 
v = |, the splitting between the two modes at the cyclotron frequency satisfies, 

at leading order in \Q\, Aui 2 = lo\ — lo 2 _ = oj 2 ^ — ui 2 , which is proportional to 
\Q\ 2 . Up to this order, experimentally one should observe one mode dispersing 
as cj+ (Eq. [|5]), and the other as ui = ui c - For v = ^ the splitting is also 
proportional to |Q| 2 . In this case one should observe both modes (u>+ and 
Eq. |^), but with different intensities. 

If the pair potential has the Coulomb form, the residue of ui- (Eq. [5l]) is 
smaller by a factor of \Q\ 2 ( m ~ 2 ^ than the residue of w+ (Eq. ^), and this is 
valid for all the values of the filling fraction that we have studied. The splitting 
between these two modes satisfies, at leading order in \Q\, Aw 2 = uo 2 ^ — uj 2 , 
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which is proportional to \Q\. For v different from i, the relative intensity 
between the two modes is proportional to (2B e g / Q 2 )( m ~ 3 > MV (Q) , which is 
bigger than 1 within our approximation. For v = the relative intensity is 
proportional to MV(Q). This factor can be written in terms of the magnetic 
length and the cyclotron energy as follows H£2iZ£, Since our approximation is 
only valid in the limit A >f, the numerator satisfies V(Q)/£ -^j-- The 

second term in this inequality is the Coulomb energy at the magnetic length, 
which is typically of the same order of magnitude than the cyclotron energy. 

Therefore, ^ ^ ~ 1- In other words, the relative intensity for 

v = ^ is also bigger than one. 



3-4 Saturation of the f-sum rule 

We show now that the long wavelength form of Koo, found at this semiclassical 
level, saturates the /-sum rule. This result implies that the non-gaussian 
corrections do not contribute at very small momentum. In Section ^ we will 
use this result to show that the absolute value squared wave function of all the 
(incompressible) liquid states has the Laughlin form at very long distances, in 
the thermodynamic limit. 

The retarded density and current correlation functions of this theory are, 
by definition 

D%{x,y) = -i9(x -y )x < G|[J„(x), J u (y)]\G > (52) 
where J M (/i = 0,1,2) are the conserved currents of the theory defined by 



Eq. 13, and \G > is the ground state of the system. Using this definition 
and the commutation relations between the currents, one can derive the /- 
sum rule for the retarded density correlation function Dq . In units in which 
e = c = h = 1, it states that 



The polarization tensor Koo and the density correlation function Doo are sim- 
ply related by Kqo — —-Doo- From Eq. [37] and [38| we see that the leading 
order term in Q 2 of the zero-zero component of the electromagnetic response 
is given by 

^0" = -T7 2 2 , ■ ( 54 ) 
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where we have used that -S = Thus, we see that the leading order term of 
K 00 saturates the /-sum rule, Eq. |5^. 

It is important to remark that the coefficient of the leading order term of 
Kqo can not be renormalized by higher order terms in the gradient expansion, 
nor in the semiclassical expansion. In the case of the gradient expansion, it is 
clear that higher order terms have higher order powers of Q 2 , and then, do not 
modify the leading order term. In the case of the corrections to Kqo originating 
in higher order terms in the semiclassical expansion, they also come with higher 
order powers of Q 2 . The reason of that is essentially the gauge invariance of 
the system. This implies that the higher order correlation functions must be 
transverse in real space, or equivalently they have higher order powers of Q 2 
in momentum space. Being higher order terms in the Q 2 expansion they can 
not change the leading order term. 

As we have already mentioned, these results hold for any model Hamilto- 
nian for the two-dimensional electron gas (2DEG) with reasonably local inter- 
actions, i. e., with pair interactions that obey Q 2 V(Q) — > as Q 2 — » 0. 



3.5 Hall Conductance, Effective Low Energy Action and Quantum Numbers 

We show now that, already within our approximations, this state does exhibit 
the Fractional Hall Effect with the exact value of the Hall conductance. In or- 
der to do so, we will calculate the Hall conductance of the whole system. Since 
we are only interested in the leading long-distance behavior, it is sufficient to 
keep only those terms in Eq. ^ which have the smallest number of derivatives, 
or in momentum space, the smallest number of powers of Q. Therefore, from 
Eq. |36| and [[t], we see that the leading long distance behavior (i. e., small 
momentum) of the effective action for the electromagnetic field is governed by 
the Chern-Simons term. In this limit Eq. ^ turns out to be 

S^(A^ « -- J ^^(-w, -Q)K x {u, Q)e^ x Q x A u (cu, Q) (55) 

where Q° = u> and Q l = — Qj. 

To study the Hall response of the system, we will now consider the limit 
of small u> and small Q. We have checked that in this theory the two limits 
commute 

Urn limi^, Q) = n ;f; 0) (56) 
lira lim^^QV ni ]°' 0) (57) 

U ^Q^0 + 111 (0,0) 
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This is a consequence of the incompressibility of the ground state. Since 9 = 
1 

2-ir 2s 



±- and Hi (0,0) = £ 



K x -> ^ = cff (58) 

Z7T 

where v is the filling fraction. The electromagnetic current J M induced in the 
system is obtained by differentiating the effective action S c s(A t _ l ) with respect 
to the electromagnetic vector potential. The current is J M = ^-e^ v \F uX . 
Thus, if a weak external electric field Ej is applied, the induced current is 
Jk = QcS^ikEi. We can then identify the coefficient 9 c q with the actual Hall 
conductance of the system a xy and get 

cr xy = OeS = 7T ( 59 ) 

which is a fractional multiple of ^- (in units in which e = fi = 1). Thus, the 
uniform states exhibit a Fractional Quantum Hall effect with the correct value 
of the Hall conductance. 

We can also derive these results by considering an effective action for the 
low energy excitations. At the gaussian level, the effective action is obtained 
by keeping only those terms in Eq. which have the smallest number of 
derivatives and we can then approximate S e s by 

cr° 
S eS ™^ L S CS (a lJt ) + -Scs(^-\). (60) 

where the induced Hall conductance er° , = ±#-. Although we have obtained 
this result within the gaussian approximation, it should be stressed that a xy 
does not get corrected to any order in perturbation theory. This is so because, 
at least for a system with a gap and,^ the thermodynamic limit, a® y is de- 



termined by a topological invariant ES^Ea, the First Chern Character C. In the 
problem at hand, we have C — p. Thus, the Hall conductance that we will 
calculate below is actually exact. 

Upon integrating over the statistical vector potentials, we find that the 
effective action for the electromagnetic fields is 

S^(A,) « 9 -fs cs (A,). (61) 

where 9 c g is given by 

ir = l + ^r- (62) 

i. c, the Chern-Simons coupling constants are added "in parallel". 
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The values of 6 and a% y determined above yield the result 

8 cS = ^ (63) 

where v is the filling fraction. As before, we obtain the value of the Hall 
conductance from the electromagnetic current induced in the system. The 
result is a xv = 8 c (r — Notice that the coefficient a°^, of the effective action 
represents the integer Hall effect of the bound states and it is different from 
the true Hall conductance a xy . 

Finally we discuss the charge and statistics of the topological excitations, 
the quasiparticles and quasiholes. We recall that, at the level of the Aver- 
age Field Approximation, the single particle excitations are fermions, single 
(charge one) composite fermions. However, the fluctuations of the statistical 
gauge fields modify this picture in an essential way. For instance, if we want to 
compute the effective charge of these excitations, we need to find the exact low 
energy behavior of the composite fermion propagator. Since, even at the level 
of the AFA, these single particle states have a finite energy gap, the Feynman 
path integral form of their propagator is dominated by smooth paths T. The 
fermions couple to both, the statistical and electromagnetic fields as if they 
had charge one. It was shown in referenceEd that the leading effect of the fluc- 
tuations of the gauge fields can be found by computing the expectation value 
of a Wilson loop operator using the effective action of Eq. |6(| More explicitly, 
we find 

,ij r <te (l (a"( a 0+A"( a 0)v = e i 7 (r) ,ij v dx„ a*(VK /g^ 
Here 7(r) is the Aharonov-Bohm phase factor 

where $o(r) = BL(T) 2 is the flux through the area enclosed by V and <f>o is the 
flux quantum. On the l.h.s. of this equation a M is the complete Chern-Simons 
gauge field, i.e., the sum of the solution of the SPE plus the fluctuating part 
a M . Eq. |6^ tells us that the effective charge of this excitation is q c g = 2s p±i ■ 
The additional phase factor can, in principle, give rise to a fractional spin but 
it is not possible to determine its value (or even if it exists!) without making 
contact with the microscopic theory. 

The statistics of these excitations can be computed using a similar type of 
reasoning. Now we need to look at a two-particle propagator and, consequently, 
we need to compute expectation values of two Wilson loops. In section |^ we 
gave a prescription for the computation of the expectation values that we need 
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here. The effective value of 9 in Eq. [Tl] now is 9 = 9 + a® y . Thus, the statistical 
angle S (measured from fermions) is 

5 = i = ^— ^- (66) 
26 2sp±l y ' 

Hence, at long distances and at low energies, the composite fermions are ac- 
tually fractionally charged anyons. This is the well know result of Arovas, 
Scrhicffer and WilczekEJ. 



4 Fermionic Chern-Simons theory for the Fractional Quantum Hall 
Effect in bilayers and partially polarized systems 

If one allows for the presence of new degrees of freedom, a richer variety of states 
can be found. The two obvious possibilities that one can consider are systems 
in which the electronic spin is not frozen by the Zeeman energy, and systems 
in which two or more layers of 2DES .axe coupled together. For instance, the 
experimentally obsepied v — I state E3, has been explained theoretically by 
Haldane and RezayiEj using the fact that the system is not spin-polarized. 

Due to continuing advances in material-growth techniques, it has been pos- 
sible to fabricate high-quality multiple 2DE layers in close proximity. In these 
systems the layer index is the new degree of freedom, and the interplay between 
the intralayer and the interlayer Coulomb interactions gives, rise to very inter- 
esting physics. In particular, this competition can explained the experimental 
observation^ of the destruction or weakening of the IQHE at odd filling frac- 
tions. Another interesting case is the one of the v = \ state. In single- layer 
systems, even though many transport anomalies have been reported—there is 
no evidence of FQHE. On the other hand, this is a well observed E3 FQHE 
state in double-layer systems. 

There are two energy scales that play a very important role in bilayers. 
One is the potential energy between the electrons in different layers, and the 
other one is the tunneling amplitude between layers. In the case in which 
tunneling between the layers may be neglected, the number of particles in 
each layer is conserved. The system has an effective t/(l) ® U(l) symmetry 
reflecting the separate conservation of charge in each layer. Tunneling processes 
break this symmetry down to a global U(l) invariance. In the next subsection 
we will describe the U(l) <S> U(l) fermionic Chern-Simons theory for bilayers 
which ignores interlayer tunneling effects. In a subsequent section we present 
an SU(2) (g> U(l) fermionic Chern-Simons theory. This theory describes the 
physics of spin as well as tunneling effects in bilayers. 
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4-1 Fermionic Chern-Simons theory for double layer FQHE systems 

The action for a double-layer 2DES in the presence of an external uniform 
magnetic field B perpendicular to it is given by 

~ J d 3 z f d 3 z> ]T W a (z)\ 2 -p a )v afi (\?-?\){\MJ)\ 2 -Pf>) 

a, (3 

(67) 

where the indices a = 1, 2 and f3 = 1, 2 label the layers, p Q is the average 
particle density in the layer a, ip(z) a is a second quantized Fermi field, /i Q is 
the chemical potential and is the covariant derivative which couples the 
fermions to the external electromagnetic field A^. In what follows we will 
assume that the pair potential has either the Coulomb form, i. e., 

V(\r\) afj = q l (68) 



r' + d (1-M) 

(with d the interlayer separation), or that it represents a short range interac- 
tion such that in momentum space it satisfies that V a p(Q)Q 2 vanishes at zero 
momentum. This includes the case of ultralocal potentials (i e., with a range 
smaller or of the same order as the cyclotron length £), in which case we can 
set V(0) — 0, or short range potentials with a range longer than £ such as a 
Yukawa interaction. _ 
Using a generalization of the methods of Section g, it is possible to showED 
that this system is equivalent to a system of interacting electrons coupled to 
an additional statistical vector potential (fi = 0,1,2) whose dynamics is 
governed by the Chern-Simons action 

S cs =Y.^fj d ** ^ VX WA (69) 

a/3 

provided that the CS coupling constant matrix satisfies 
a g _ 1 ( 2s 2 -n 



K ° ~ 2tt(4 SiS2 - n 2 ) { -n 2s x ) (70) 

where si, S2 and n are arbitrary integers. In Eq. |6^ xq, x\ and x 2 represent the 
time and the space coordinates of the electrons respectively. In the equivalent 
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theory the covariant derivative is given by 

D% = d li + i~A ll + ia% (71) 

and it couples the fermions to the statistical gauge fields (a^), and to the 
external electromagnetic field (A^). Notice that the theory has now a U(l) ® 
U(l) gauge invariance. 

The Chern-Simons coupling constant matrix k"' 5 of Eq. [70] must be non- 
singular. For the case of a singular matrix, it is better to arrange the gauge 
fields in linear combinations a:£ which diagonalize the coupling matrix. In 
general this situation arises only if the matrix elements satisfy Kn = K22 = k 
and K12 — K21 — k' and the matrix is singular if k = n' . In this case the 
gauge field at acts like a local constraint on the allowed states demanding 
that pi(x,t) = p2(x,t). In other words, the allowed states consist of pairs 
of fermions, one from each layer. If we think of the layer index as spin, these 
states are local spin singlets and the allowed states are spinless bosons of charge 
2 (with zero transversal spatial extent). Thus, for a singular coupling matrix, 
this theory effectively describes a FQHE of a single layer of charge 2 bosons. 
While this is prefectly correct, it is clear that this is an extreme case and that 
this approach is uncapable of describing a singlet FQH state not made of local 
singlets. To describe a state of this sort within a Chern-Simons approach it 
is required to attach separate fluxes to each species of fermions in a manner 
consistent to the full SU(2) symmetry of spin. To describe this situation 
correctly it is necessary to use a non-abelian approach, which we present in 



Section 4.4 



There is another limit of physical interest. Clearly, in a bilayer system 
(or for a single layer with spin) it should also be possible to attach the same 
flux (using only one gauge field) to the fermions of both layers. This situation 
is achieved by taking the limit k — n' — > 00 (which kills the field at) while 
holding + k!) finite. This_is the conventional scheme used to describe 
partially polarized FQH states E3. 

The Chern-Simons action implies a constraint for the particle density jg (x) 
and the statistical flux B a = e^c^a", given by 

]S{x) + KaP B p {x) = Q (72) 

(from now on we assume that repeated indices are contracted). 

This relation states that the electrons in plane a coupled to statistical 
gauge fields with Chern-Simons coupling constant given by Eq. [7^, see a sta- 
tistical flux per particle of 27r2s a for the particles in their own plane , and 
a statistical flux per particle of 2im for the particles in the opposite plane. 
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(Notice that in units in which e = c = K = 1, the flux quantum is equal to 2tt) . 
If the coefficient of the Chern-Simons term is chosen with the above prescrip- 
tion, all the physical amplitudes calculated in this theory are identical to the 
amplitudes calculated in the standard theory, in which the Chern-Simons field 
is absent. Of course, this is true provided that the dynamics of the statistical 
gauge fields is fully taken into account exactly. 

Here we take into account the dynamics of the Chern-Simons gauge fields 
in a semiclassical expansion, which is a sequence of well controlled approxima- 
tions. In practice, we will only consider the leading and next-to-leading order 
in the semiclassical approximation. 

Using the constraint enforced by the Chern-Simons action, the interaction 
term of the action Eq. |67| becomes 

Sint = -\J d3z / dV ^ aS Bs(z)~p a ) V a p(\z-z'\) {k^B 7 {z') - p p ) (73) 

The quantum partition function for this problem is, at zero temperature 

Z[AJ = I Vrv^Va" exp(iS(r^,a%A^) (74) 



Since the action is quadratic in the fermions, they can be integrated out. The 
effective action (<S c ff) is given by (in units in which e = c = h = 1) 




(75) 

where 

- A«) = ~J d 3 z J d 3 z' (K a \B s (z) - B 5 (z)) p a ) 

V a0 (\z- zi\)(K^(B 7 (z') - B^z')) - pp) (76) 

Here we have written the external electromagnetic field as a sum of two terms, 
one representing the uniform magnetic field B, and a small fluctuating term 
A^ whose average vanishes everywhere. The latter will be used to probe the 
electromagnetic response of the system. Notice that we have used the in- 
variance of the measure £>a" with respects to shifts, to move A" out of the 
covariant derivatives and into the Chern-Simons and the interaction terms of 
the effective action (Eq. ffsl ). 

The path integral Z can be approximated by expanding its degrees of 
freedom in powers of the fluctuations, around stationary configurations oiS e ff. 
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This requirement yields the SPE or classical equations of motion. 

<jS(z)>f = - KaP [<B f3 (z)> -B?{z)>] 
<3%(z)> F = -K a peia[<ef(z)>-E?(z)>] 

+e k idf J d 3 z' K a 'V eS (z, z')[-k 5 ^ <B"i-B~i> (z r ) - p 5 } 

(77) 

Just as in the case of the single layer problem, these equations have many 
possible solutions. The uniform solutions of Eq. |77] represent the FQH fluid 
states of the bilayer system at the mean field level. In this state, the electrons in 
the layer a see a total flux B" s , equal to the external magnetic flux partially 
screened by the average Chern-Simons flux, i. e. B" s = B + (B a ) = B — 
(k~ 1 ) q/3 P / 3. It is easy to see that the uniform saddle-point state has a gap only 
if the effective field B" s is such that the fermions in layer a fill exactly an 
integer number p a of the effective Landau levels, i. e., those defined by B^ s . 
In other words, the AFA to this theory yields a state with an energy gap if the 
filling fractions of each layer satisfy 



f 



v\ = 



Pi 1 Ni 



V2 = ±± + 2s 2 + n fr (78) 

where N\ and N 2 are the number of particles in layers 1 and 2 respectively. The 
sign in front of p\ and p 2 indicates if the effective field is parallel or antiparallel 
to the external magnetic field. 

Using the fact that v = v\+v 2 , and that the number of flux quanta enclosed 
by each plane is the same, the filling fractions can be written as follows 



(±^ + 2*2) 



-(±± + 2 Sl )(±j- 2 +2s 2 ) 



" 2 = ni-(±± + 2 Sl )(±j- 2+ 2s 2 ) (79) 

where p\, p 2 , s\, s 2 and n are arbitrary integers. For the special case in which 

the two layers have the same occupancy, N\ = N 2 and i>\ — v 2 = the 
allowed fractions are 

2p 

v{p, n, s) = (80) 
[n + is)p + 1 
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where p is an arbitrary (positive or negative) integer. 

The effective magnetic field can be written in terms of the external mag- 
netic field as 

% = (81) 

For general values of these integers, the states whose filling fractions are given 
by Eq. (79] have a gap at the mean field level of this approximation. This en- 
sures that the perturbative expansion is meaningful. If there is no gap for the 
excitations of the mean field ground state, the perturbative expansion breaks 
down. The breakdown is signalled by infrared divergencies at low tempera- 
tures. Such is the case for the compressible or "Fermi Liquid" states. 



It is clear from Eq. 79 that in the bilayer systems, there are many possible 
choices of the numbers pi, p 2 , s±, S2 and n, i. e., many different states, which 
have the same filling fractions. These states are physically different since, as 
we will see, the quantum numbers of their excitation spectra are different. This 
is an example^ of the concept of topological order that has been formulated by 
X. G. Wen l3. More recently, HaldaneEll has proposed an algebraic criterion 
to determine what states are topologically stable. We will see that all the 
solutions with an energy gap for all excitations are topologically stable in the 
sense of Haldane. 

As a result of the presence of many solutions, the phase diagram for bi- 
layer s is much richer than for spin polarized electrons in a single layer. Which 
particular solution is actually realized is determined by the microscopic dynam- 
ics of a particular system and it requires a detailed calculation of the ground 
state energy beyond the^verage Field Approximation. Experiments on spin 
polarized 2DEG bilayersE3c3 at filling fractions v = 1, 1/2 have shown a non- 
trivial phase diagram with at least two phases: a compressible phase in the 
regime in which the 2DEG's are well separated and an incompressible_uhase 
when they are closer by. However, as it was emphasized by Wen and ZeeEil, the 
incompressible states at these two filling fractions actually have quite different 
properties. We will see now that the phase diagram can be quite complex. 

Let us consider some cases of special interest, in particular with Ni = N2- 
For example, at level one of the hierarchy, we choose p\ = p% = 1,. flfijjr-J = m\ 

J LIU -\ 



and 2s2 + 1 = UI2, and we obtain the so called (mi, mj, n) statesEi'tSea whose 
filling fractions are v — "^^^^2" ■ I n particular, for mi = to 2 = 3 and n = 1, 
this is the (3, 3, 1) state, whose filling fraction is v = ~. For-n^^p p 2 = 1, and 
2si + 1 = 2si + 1 = n = to, we obtain the (m, to, to) states whose filling 
fractions are v = — . 

The states (mi, TO2, n) are the equivalent of the Laughlin states for bilayers 
and as such should be regarded as primary states. Hence, for each (mi, TO2, n) 
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primary state there is a hierarchy of states labelled by the index p. For example, 
the primary state (1,1,1), defined by p — 1, s = and n = 1, has the hierarchy 
v+{p, 1, 0) = = 1, |, |, |, . . p. 2. The first three states in this series have 
already been seen in experimentsCJ. Similarly, the primary state (3, 3, 1) (which 
has p = 1, n = s = 1), has the hierarchy v+{p, 1, 1) = g^px = 5, f , f , ■ ■ ■ , § • 
Recall that v±{p 1 n, s) is the io£a/ filling fraction of the bilayer and that the 
individual filling fractions of each layer are v/2. 

We can also consider the limit p\ = P2 = p — » 00. In this case, if s± — S2 — 
s, the state with filling fraction v — i will be obtained for all the values of n 
and s such that n + 2s = 2m. Since p — > 00, the effective field vanishes, and 
we find the analogous of the^ompressible states for the single layer problem 
discussed by Halperin et al Ex All of these states are degenerate with the 
v = — states mentioned above. 

rn 

Clearly, it is always possible to construct a large family of, in principle, 
distinct states which have the same filling fraction. In particular, new states 
are generated by transferring mira-plane flux-particle attachment, determined 
by si and S2, to mter-plane flux attachment, determined by n. Since the inter- 
plane attached fluxes are negative in these states, the particles of one plane see 
the particles of the other plane as if they were holes. Thus, in these states, there 
is an effective attractive force between particles on different planes. Therefore 
the wave functions, instead of having zeros of higher order (which represent 
repulsion), should have a larger weight when particles from different planes 
are closer to each other (in the sense of their coordinates projected on the xy 
plane). Asymptotically, these wave functions appear to have "poles" in the 
interlayer coordinates. One could think of this state as a condensate of inter- 
layer bound states. Such states should become preferable as the interlayer 
coupling constant becomes large. 



.2 Spectrum of bilayer FQH states 



Following the same steps as in Section 3.1, we can integrate out the gaussian 
fluctuations of the statistical gauge fields a", and obtain the effective action 
for the electromagnetic fluctuations A" t 

S*W) = ljd 3 xj d 3 yA^x)K^(x,y)A^y) (82) 

Here K^g is the electromagnetic polarization tensor, which measures the lin- 
ear response of the system to a weak electromagnetic perturbation. As for 
the single-layer systems, we can use this effective action to calculate the full 
electromagnetic response functions at the gaussian level. 
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In this Section we discuss the spectrum of excitations for two different 
states, the (to, to, n) and the (to, to, to) states. The derivation of these results 
can be found in reference Ell. 

We first discuss the case of the (to, to, n) states. For simplicity, we study 
the collective modes for the state (3,3, 1). All the other states can be studied 
by straightforward application of the same methodsEH. 

In this case the total filling fraction is v = i, being v\ = v 2 = \ ■ The 
effective cyclotron frequencies and magnetic fields are u)^ s = = tl> = ^ 
and Bl s = Bl s = B = f . 

We find that there is a family of collective modes whose zero-momentum 
gap is kuj, where k is an integer number different from 1. At mean field level, 
there are two modes for each integer multiple of ui. After including the gaussian 
fluctuations we find that there are no modes with a zero momentum gap at 
u>. One of them has been "pushed up" to the cyclotron frequency and the 
other up to 2ui (at Q = 0). Therefore, at these multiples of u> there are three 
degenerated modes for Q = 0. For Q ^ 0, the degeneracy is lifted and these 
three modes have different dispersion curves. 

At 2<D we find that there are two modes with residue Q 2 , and one with 
residue Q A . The former are 



i())= (2^) 2 + (-|^)^ 2 a ± 



(83) 



where 



3M . 

Z7T 



v 12 ) 2 + 16 



(84) 



Here v a p are the zeroth order coefficients of the Fourier transform of the inter- 
particle pair potential for short range interactions. For Coulomb interactions 



and v 12 



r e -\Q\d 

6 ra8 



if \Q\d <C 1, therefore a± = ±4 in this limit. 



The residues in K q corresponding to these poles are 



Res(K$,u ± (Q)) 



0J C 



8tt 



1 



1 



-1 



1 



(85) 



It is clear from the form of these residues, that these excitations are out of 
phase modes, because they only couple the "out of phase" density (i. e., they 
couple p~ with itself). 

The other mode at 2u> is 



zi>cff 



(86) 
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and its residue is proportional to Q 4 ( } \ I . Thus, this is an in phase 



1 1 

mode since it only couples p + with itself. 

There are three modes whose zero momentum frequency is the cyclotron 
frequency. The only mode with residue Q 2 has energy 

= uj c + (2 H )3( jo; (87) 



and residue 



Res(K 00 ,oj(Q)) = -Q 2 w c £- (J J), (88) 

As in the case of single layers, except for the lowest branch of collective modes 
and the cyclotron mode at Q = 0, all the other modes become damped due 
to processes induced by the non-quadratic interactions among the collective 
modes. 

We now discuss the collective modes of the (to, to, to) states. In this case 
the total filling fraction is v = ^, being v\ = v 2 = ^ • The effective 
cyclotron frequencies and magnetic fields are = lu 2 s = uj = ^ and 

d1 n2 — d S 

-^eff - B eff = B eS ~ 2^- 

We find again a family of collective modes whose zero-momentum gap is 
ku), where k is an integer number different from 1. At mean field level, there 
are two modes for each integer multiple of u>. After including the gaussian 
fluctuations we find that there are no modes with a zero momentum gap at u>. 
One of them has been "pushed up" to the cyclotron frequency. Therefore, at 
u) c there are three degenerate modes for Q = 0. For Q ^ 0, the degeneracy is 
lifted and these three modes have different dispersion curves. The other mode 
at u has been "pulled down" to zero frequency at Q = 0, i. e., it has become 
a gapless mode. 

We will distinguish between the cases to = 1 and m/1. 
a) Case to = 1 

The gapless mode is 

lu{Q) = v s \Q\ (89) 

where 

M 

where v a p are the zeroth order coefficient of the Fourier transform of the in- 
terparticle pair potential for short range interactions. For Coulomb interac- 
tions Vn{Q) = ^ and v 12 (Q) = ^- e ~^ d w 4 if \Q\ d < L Thcrcforc 
(un — U12) = for Coulomb interactions (in the limit \Q\d <C 1). 



l + ^(«u -^12)]^ (90) 
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The residue in Kqq corresponding to this pole is 



Res(K 00 ,u(Q)) = -Q 2 u c 



1 -1 
-1 1 



(91) 



Therefore, this is an out of phase mode. 

At lo c = 2u) we find that there are two (in phase) modes with residue Q 2 



W±(Q)= uj 2 c + (—— Q 2 a± 

L Z-D e ff 



where 



M, , ,!.„ .„ 

OL± = 2^(^11 +Vi 2 )±\l (— ) + V12Y + 16 



(92) 



(93) 



For Coulomb interactions t>n + V\2 = 2^- if \Q\d <C 1, therefore this term is 
higher order in Q and it should be neglected, i. e., a± = 4. 



The residues in Kqq corresponding to these poles are 



Res(K 00 ,uj±(Q)) = -Q 2 iu c 



8- 



1 1 
1 1 



16 



therefore these are in phase modes. The other mode at lu c is 

^o(Q) = 



;2 _ 2 (-^-) CO 2 



(94) 



(95) 



and its residue is proportional to Q 4 



(out of phase mode). 



1 -1 

x -1 1 
b) Case m/1 

The gapless mode is an out of phase mode with the same form as for m = 1 
(Eq. ^) and with the same residue (Eq. 

At uj c = 2mu) we find that there is only one in phase mode with residue 

Q 2 



>(Q) = 



2m - 1 M 



771 — 1 



2tt 



Q 2 



'2B, 



err 



The residue in Kqq corresponding to this pole is 



Res(K Q0 ,Lu(Q)) = -Q 2 uj c 



8n 



1 1 
1 1 



(96) 



(97) 



The modes with zero jnomentum frequency u> — kuj for k > 3 have been 
calculated in reference E3. 

All the considerations about the validity of this spectrum of collective ex- 
citations beyond the semiclassical approximation that we discussed in previous 
sections are of course valid in this case as well. 
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4-3 Effective Low Energy Action, Hall Conductance and Quantum Numbers 



We will derive now the effective action for the low energy degrees of freedom 



of the bilayer system. Using the same line of argument of Section 3.5, we keep 
only the leading terms in a gradient expansion and find 



S cS (5« A°) « ^- / ' o?x a% e^ vX d x ~4 - ^- / d 3 x a« e^ vX d x A? 

~ J d 3 x A° e^ x d x ~4 + ^ J d*x A° e^ x d x A? 



(98) 

where we have set k a/3 = %r-5 al3 + n a ^ . 

The effective action of Eq. B$ is very similar to the effective action of Wen 
and Zee's Matrix classification of Abelian quantum Hall states c3. It differs 
from it in two respects: in Eq. |98| there is an induced Chern-Simons term 
for the background electromagnetic fields, and the coefficients are not quite 



the same as in reference CJ. However, as we will see below, Eq. 98 leads to the 
correct value of the Hall conductance and of the fractional charge and statistics 
of the quasiparticles. Thus, they describe the same physics. 

The next step is to integrate the statistical gauge fields to obtain the 
effective action for the electromagnetic field. In particular, we will need to 
compute the inverse of the matrix R a ^ . This inverse only exists if A = [(^- + 
2si)(^- + 2s2) — n 2 ] ^ 0. Therefore we must consider two cases, the one in 
which A^O and the one in which A = 0. 
Case A ^ 

Upon integrating over the statistical gauge fields in Eq. ^8[ the effective 
action for the electromagnetic field results 

St$(A«) « ^ J d 3 xA« e,„ x d x ~Al (99) 

where the effective quantum Hall conductance matrix 

k££ = K Q T[^' 3 - (k-^V 3 ] (100) 

is given by 

K «fi 1 ( (£ + 2 S2 ) -„ \ 

~ 2*[(£ + 2i)(£ + 2s 2 ) **] { -n (i + 2«x) ) [W1) 
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The matrix gives the intralayer and translayer conductances of the bilayer 
system. The translayer conductances characterize the Coulomb drag effects in 
bilayers. 

In particular, let us consider the case in which both layers are coupled to 
the same electromagnetic field. Then = A^ = A^, and the coefficient in 
the effective action results 

El 2n - (±- + 2si) - (— + 2s 2 ) v 
R cS ~2tt n ^(±- + 2 Sl ){±- + 2s 2 ) ~ 2lr ( 1U2 J 

a/3 v Pl L,y P2 Al 



Following the same arguments as in Section |3.q we find that the electromagnetic 
current induced in the system is J M = ^ L e f2iy \F^ x . Thus, if a weak external 
electric field Ej is applied, the induced current is Jk = K e seu~Ei. Therefore 
the coefficient k c r is the actual Hall conductance of the system. 



Vxy = Kcff = 7T~ (103) 
Z7T 

which is a fractional multiple of ^- (in units in which e = Ti = 1). Thus, the 
uniform states exhibit a Fractional Quantum Hall effect with the correct value 
of the Hall conductance. 

The effective action of Eq. |9^ also determines the quantum numbers {i. e., 
charge and statistics) of the vortices (topological excitations) supported by 
this FQH state. These quantum numbers were calculated in reference E3 and 
we will not give an explicit derivation here. The charge of the excitations was 
found to be given by 

<& = £ (104) 

In particular, for the (m, m, n) states, the effective charge in both layers is the 
same and is given by q cS = 

On the other hand, the statistics of the excitations on the same layer is 

S a =n(l + = fo±gfej a ) (105 ) 

" V /„ 2s 2 \f„ i 2si \ n 2 I v ' 

\Pl -T 4 SlS2 - n 2)(P2 -T 4siS2 _„2j (4sis 2 -n 2 ) 2 

where if a = 1, then [3 = 2 and viceversa. In particular, for the (m,m,n) 
states, the statistics is 8 = — ^Fl^ , independent of the layer. The relative 
statistics (i. e., for excitations on different layers) 612 is given by 



5x2 = TT ^ iSlS2 ~f Sl -2 (106) 
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In particular, for the (to, to, n) states, the relative statistics results 5 = r^rz^s% ■ 
It is worth to-discuss these results in the context of Haldanes's topological 
stability criterionc3. After some (tedious) algebra it is straightforward to show 
that Haldanes's topological stability is equivalent to the requirement that the 
denominators in Eq. 105 be non-zero. For the special case of the (to, to, n) 
states, this holds if to ^ n. We will see next that this requirement fails if there 
exists one branch of collective excitations with a gapless mode, such as in the 
(to, to, to) states. 
Case A = 

It can be shown that when R a P is not invertible, i. e., it has a zero eigen- 
value, the corresponding linear combination of the gauge fields become mass- 
less. In other words, the in phase gauge field a+ — a* + ai has a finite gap 
which couples to the electromagnetic field A^, while the out of phase gauge 
field a~ — — ai is gapless. 

We will study in particular the case of the (to, to, to) states which satisfy 
the condition A = 0. For these states 



1 1 

2tt(1 - 2m) V 1 1 



(107) 



which is clearly non invertible. 

We can write the effective action for the fluctuations of the Chern-Simons 
gauge fields (Eq. p8j) in the basis defined by = a* ± a 2 



SM~a«,A«) « -\ J ^ at 

~\ J d 3 x — - — {-~a+e^ x d x A+ + (2m - l)a-e^ uX d x A~) 

~2 J 4tt(1 - 2m) (~^ c ^^ ~ a - + ( 2m ~ ^)K^uxd x a~) 
+ \ J d 3 x 4?r(1 [_ 2m) {-A+e^ x d x At + (2m - l)A~e^ A~) 

(108) 

The gauge field a~ appears as a Lagrange multiplier in this action. The in- 
tegration over it states that the current e^, v \d x A~ vanishes. This is trivially 
valid if the electromagnetic field is the same for both layers, because A~ = 0. 

The integration over at gives the effective action for the field At- If we 
consider the case in which A}. = Af, = A u , the result is 

S e T = Jd 3 xA^ vX d x A u (109) 
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Following the same steps as in the previous case, the Hall conductance results 
° x v = l^m = 5tF which is the correct value for the (m, to, to) states. 

The existence of a gapless branch in the spectrum of collective excitations 
has important consequences for the spectrum of topological excitations i. e., 
quasiparticles and quasiholes. It was shown in reference E3 that the gapless 
collective excitations mediate an effective confining gauge force between topo- 
logical excitations which carry the quantum numbers of the out-of-phase (or 
inter-layer) mode. The result is that the only surviving excitations are bound 
states which are neutral with respect to the inter-layer gauge group but carry 
charge 1/m and fractional statistics ir/m (measured from fermions). These 
states have a finite energy gap. For the particular case of the state (1, 1, 1), 
these bound states are bosons of charge 1 and finite gap. This gapped bound 
state has the same quantum numbers as the skyrmion excitation of the (1,1,1) 
discussed, by K. Yang et. al, who view the (1,1, 1) state as a quantum ferro- 
magnetic. S. L. Sondhi and collaborators had predicted earlier that a partially 
polarized v = 1 QH state should have skyrmion states E3. Our results indi- 
cate that all of the (to, to, to) states should have skyrmion excitations but, in 
general, they have fractional charge X/m and fractional statistics ir/m. Hence, 
the bilayer state (3,3,3), which occurs at v = 1/3 (1/6 on each layer), has 
skyrmions of charge 1/3 and statistics tt/3. 

The reduction of the number of topological excitations for the case m = n 
illustrates, in a physically transparent way, the meaning of Haldane's crite- 
rion for topological stability. Indeed, Haldane's criterion predicts that the 
(to, to, to) states are not topologically stable. This of course does not mean 
that the state itself is unstable or unphysical. Rather, it means that if one 
or more branches of collective excitations become gapless, the quasiparticles 
which couple to them are not present in the spectrum and the corresponding 
quantum numbers are not allowed. 



4-4 Spin Singlet States and SU(2) Symmetry 

In this Section we discuss the application of Chern-Simons methods to a 2DEG 
with an SU (2) symmetry. This symmetry may_be exact or may be broken by 
explicit terms in the Hamiltonian. In referenceca a non-abelian Chern-Simons 
approach was used to construct a theory of the spin singlet (to + 1, m + 1, to) 
Halperin states. Here, we present a generalization of the approach of reference 
C3 which will yield the SU{2) hierarchies. A related, and completely equivalent 
approach was developed independently by Frohlich, Kerler, and MarchettiCJ. 

In the approach of Balatsky and FradkinHa (BF), instead of the two- 
component abelian Chern-Simons theory that we use in this Chapter, a non- 
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abelian Chern-Simons gauge field is introduced. The advantage of the BF ap- 
proach is that the SU(2) invariance is manifest and it is not the consequence 
of a subtle dynamical mechanism. The disadvantage of the BF approach is 
that the non-abelian Chern-Simons theory is substantially more sophisticated 
and technically more demanding than the abelian theory used before. In the 
BF approach, the electron is viewed as a composite object which is made of a 
particle that carries the charge (the holon) and another particle that carries 
the spin (the spinon). This arbitrary separation gives rise to the existence of 
an abelian gauge symmetry (called RVB by BF). The requirement of gauge 
invariance forces the holons and spinons to be glued together in bound states, 
the electrons. In the BF approach, the need of a non-abelian gauge field is a 
consequence of the assignment of fractional (semion) statistics to both holons 
and spinons. In this way, SU (2) fluxes are attached to a set of charge neutral, 
spin-i, fermions, which become the spinons. The SU(2) symmetry only ad- 
mits Bose, Fermi or semionic statistics. The holons, instead, are represented 
by fermions attached to U(\) fluxes. The U(l) and SU(2) Chern-Simons cou- 
pling constants must be chosen in such a way that the holons and spinons 
are semions. Within this approach, the FQHE of the spin-singlet states is the 
FQHE of the semions. The spin structure just sits on top of the FQHE. 

There is another, and more obvious, way to attach fluxes to particles while 
keeping the SU{2) spin symmetry untouched. Belkhir and JainE2l proposed 
a spin singlet wave function for the state with filling fraction v = i-, based 
on a composite fermion picture. In their construction, they attach the same 
number of pairs of flux quanta to both up and down spins. These approaches are 
represented in the U (1) ®U(X) theory of the preceeding sections by demanding 
that both up and down electrons see that same flux at all times. This means 
to choose a Chern-Simons matrix K a/3 which is proportional to the identity, 
i. e., s\ — S2 and n = 0. It is easy to use this approach to get the SU{2) limit 
of the Jain states, but with arbitrary polarization. However, it is not a very 
efficient approach to get all of the spin singlet states. The approach of BF deals 
with this states more directly. Alternatively, it is possible to attach flux to the 
charge degrees of freedom, leaving spin untouched. This construction, but in 
the bosonic Chern-Simons language, was used by D. H. Lee and C_KaneEil 
and subsequently applied by Sondhi, Karlhede, Kivelson and RezayiEl in their 
theory of skyrmion states in polarized FQHE states. This approach does not 
describe the spin singlet states. In what follows, we follow the BF approach. 

More concretely, we introduce a holon field (f> and a spinon field Xa {pt =T 
, J.) and represent the electron field operator ip a as ip a (x) = </>(%) Xa( x )- BF 
showed that, for an SU(2) invariant system, the system defined by the action 
of Eq. ^3?] is equivalent to the following theory of (interacting) spinons and 
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holons 

S ^charge ~l~ ^spin ~l~ ^interaction (H^) 

where the action for the charge degrees of freedom is 

Charge - J d 3 x Ut(x) (iD c + fj,) <j>{x) + ^<t>\x) D 2 c <f>(x) 

d 3 x -e IJtvX ai 1 {x)d v a x {x) (111) 
while the action for spin is 

S spin = / d 3 x (xl{x) iD s X *(x) + 2j{Xl(x) D 2 sXa (x) 

- J d 3 x^e^MV l (xWb\x) + h^(x)b l '(x)b x (x)) (112) 

47T 3 

and an (instantaneous) pair interaction term 

interaction = - / <?X J d 3 x' ^(p(x) - p)V (x - x')( P { X ') - p) (113) 

where p(x) is the charge density 

p(x) = iPt(x)*t> a (x) ee tf[x)<t>[x) (114) 

(see below). In Eq. |lll[ a M is the statistical vector potential which turns 
the holons into semions. This condition requires that the U(l) Chern-Simons 
coupling constant 9 be restricted to the values 

1 2tt 

- = + 2tt2s (115) 

m 

where semion st atist ics requires that m = ±2 and s is an arbitrary integer. 
Likewise, in Eq. 112 is the SU(2) non-abelian statistical gauge field which 
takes values on the SU{2) algebra. Hence, we can expand the field in the form 
bfj,(x) — b^(x) r a , where r a (a — 1,2,3) are the three generators of 5/7(2) in 
the spinor representation , i. e., the set of 2 x 2 Pauli matrices. The SU(2) 
Chern-Simons coupling constant fc, the level of the Chern-Simons theory is , for 
our system, equal to k = ±1e£I. Hence, we have a ("trivial") level one Chern- 
Simons theory. All the representations of a level one Chern-Simons theory 
are known to be abelian and to correspond to abelian fractional statistics of 
fermions, bosons or semionsQ. The only ambiguity left is the sign of k which is 
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the chirality (or handedness) of the scmion. There is a similar sign ambiguity 



in the coupling constant 9 in Eq. HE. Different choices of these signs lead to 
different statistics for the bound states of holons and spinons. The requirement 
that the bound state be an electron, which is a fermion, leads to the condition 
sign(fc) = sign(m). 

The U(l) and SU(2) charge and spin covariant derivatives are 

Dl = df,- z(A M + a,, + c p ) 

D; = Idp-ifa-Icn) (116) 

where is the external electromagnetic field and / is the 2x2 identity 
matrix. The gauge field c M is the U(l) "RVB" gauge field which glues spins 
and charges together. The covariant derivatives have been chosen in such a 
way that holons and spinons have opposite charge with respect to c M . Hence, 
the strong fluctuations of this field binds holons and spinons into states which 
are locally singlets under the "RVB" gauge transformations, i. e., electrons. In 
fact, all the gauge field does is to enforce the constraint that the 3-current 
of the holons equals the 3-current of the spinons, as an operator statement in 
the physical Hilbert space. This is seen clearly from the equation of motion 
generated by c M 

jRVB W = Jl^ = o tf(x)<l>(x) - xt(x)xAx) = (117) 
Sco(x) ^ 

Thus N c , the number of charges, and iV-j- and N±, the number of up and down 



spins, must obey the obvious relation N e = JV-j- + iVj_. Also, Eq. 117 tells us 
that the local particle density operator p(x) can be identified with the holon 
charge density operator, i. e., p(x) — <f>*{x)(f){x). 

In the way the theory has been set up, it may appear that this theory may 
apply only if SU (2) is an exact symmetry. However, this is not the case. Let 
us consider for instance the effects of an SU(2) symmetry breaking field of the 
form of a Zeeman term in the action 



^Zccman 



= / d 3 x B ■^ a (x)T a ^ p (x) (118) 



Using the definition tp a (x) = 4>(x) Xa( x ) and the constraint Eq. 117, we can 
write the Zeeman term in the equivalent form 



'Zccman 



d z xB-xl{x)f M {x) (119) 



45 



Hence, the Zeeman operators of Eq. 118 and Eq. 119 have exactly the same 



form. This is an exact identity in the restricted Hilbert space of states satisfying 



the constraint of Eq. 117. Likewise, spin-spin pair interaction exchange terms 
can be written exactly in the same way. Clearly, these rules also apply to any 
term involving spins even if the symmetry is not respected. Thus, a bilayer 
2DEG is a system of fermions with a layer label which can be regarded as 
a "spin" index, even though the SU(2) symmetry is actually broken. The 
interlayer tunneling term can be regarded as Zeeman term. Similarly, the 
matrix pair interaction of the electrons in the bilayer system can be represented, 
in the SU(2) formulation, in terms of a suitable combination of the SU(2) 
invariant charge pair interaction and an Ising interaction term. 

We will not go into the details of the full non-abelian theory. Instead, 
we will use this theory to determine the allowed fractions for SU(2) invariant 
states. Thus, we will just consider the AFA equations for this theory. There 
are two sets of AFA equations, one for the charge sector and one for the spin 
sector. The AFA equations for the charge sector are just the AFA equations for 
a charged interacting semion liquid with N c particles in an uniform magnetic 
field with N$ flux quanta, i. e., a FQHE of semions. A simple application of 



the methods of Section 4.1 yields the constraint for the charge density operator 



=> j (x) = -0B c (x) (120) 



5a (x) 



where jo(x) is the charge density in the liquid state, i. e., p, and B c {x) is the 
U(l) statistical charge flux. This equation, when specialized on fluid states, 
yields the allowed fractions. Similar considerations for the SU(2) gauge field 
6° (with a — 1,2,3) yield a constraint for the local spin density operator 

=> jS(x) = -i>lr a ^ a = t-B: (121) 



8b^{xY JUV ' a 2tt 

where B^ is the SU (2) spin flux. 

For fluid states, we generalize the Average Field Approximation and re- 
place these exact local operator identities by translationally invariant averages. 
This replacement may be problematic in the case of the RVB gauge field since 



its only mission is to enforce exactly the constraint Eq. 117. At the level of 
the wave functions for the allowed ground states, this constraint simply means 
that every coordinate for a charge degree of freedom has to coincide with the 
coordinate of a spin degree of freedom. 

Thus, we seek fluid states with N e particles, Nj, flux quanta and filling 
fraction v = Ne/N^. We wish to determine the filling fractions and spin for 
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which the ground state is a fluid. In addition to the usual AFA equation for 
the charge, that will give the allowed fractions, we will now get conditions for 
the spin and polarization of the allowed states. 

In the charge sector, we have a FQHE of spinless semions. In this case, the 
AFA consists of a system of N e spinless fermions filling up effective Landau 
levels, exactly as in the case of spin polarized electrons (and anyons!). The 
same line of argument that was used to derive the allowed fractions in Section 



4.1 now tells us that, from Eq. 120, the effective number of fluxes is 



** = N *-£o (122) 

and, hence, the allowed fractions, u^, for the SU(2) fluid states satisfy 

v p 2lTV 

where p is a positive integer and the + sign corresponds to a "particle-like" 
FQHE (i. e., the effective flux parallel to the external flux) while the — sign 
holds for a "hole-like" FQHE (i. e., the effective flux anti-parallel to the exter- 
nal flux) . By using the allowed values of 9 we find that the allowed fractions 
for the SU{2) fluid states are of the form 

v ± (v s-m) = — = ^ (124) 

v ' ; (2sm - l)p + m ±2+ (4s - sign(m))p v ; 

where we have specialized for the case of interest, m — 2sign(m). 



The hierarchy of FQHE states of Eq. 124 is a generalization of the states 



found by BF, which are obtained by setting p = +1 and m = +2, i. e., 
^ + (1, s; —2) = j—tt — 2, 1, f , ■ • ■, which coincide with the Halperin spin singlet 



124, where 



states. The state with v = \ is also part of the hierarchy of Eq. 
it is realized as the state with p = +2 (—2) for m = +2(— 2) and s = 1. 
Numerical studies show that, for Coulomb-like interactions, this state is not 
favored and that the (3, 3, 1) state is an accurate representation of the ground 
state. Notice that this hierarchy includes a state at filling fraction v = | . This 
state is found as the p = —10 ("hole-like"), m = —2 and s — member of the 
SU(2) hierarchy, or as the p\ — P2 = —5, Sj = s 2 — s and n + 2s = 1 member 
of the 17(1) <g> U(l) hierarchy. 



In contrast with the states found in Section 4.1 using the abelian theory 
for bilayers, the states of the SU(2) theory at each filling fraction, are arranged 
in irreducible representations (multiplets) of SU(2). The spin and polarization 



of the states is determined by Eq. 121. The only subtlety here is that, since 
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the components of the total spin do not commute with each other, one can 
only determine the total spin and total projection. For a system without a 
boundary, the choice of total spin S and total projection along an arbitrary 
polarization axis, say S z , are constants of motion which are invariant under 
local SU(2) gauge transformations (but, of course, change under global SU(2) 
rotations). There are two generic situations of physical interest: (a) spin singlet 
states (or with microscopic total spin S/N e m 0(1 /N e )) and (b) states with 
macroscopic spin, S ~ N e , i. e., ferromagnetic states. The total z-component 
of the spin polarization M = \ (N^ — N± ) obeys 

2M = N 1 -N l = ^-{BI)L 2 (125) 

where 1? is the area. It is clear that it is possible to construct all multiplets 
with spin \S\ < ^f-. In the thermodynamic limit, the spin singlet states have 
5 = and, hence, iVj — N± — N e /2. In contrast, the ferromagnetic states 
have, with an appropriate choice of the quantization axis, a non-vanishing 
extensive value of M and, hence, a non-zero value of £>J. 

The wave functions for the spin sector of the spin singlet states have to be 
determined from the states of a level one SU(2) Chern-Simons gauge theory 

^ith N e sources in the fundamental representation. It was shown by Witten 
that these wave functions are correlation functions of conformal blocks of 
a conformal field theory in two Euclidean dimensions, the SU(2) level one 
Wess-Zumino- Witten model. This fact was used by Read and Moore E3 and 
by Balatsky and Fradkinc3 to show that, the wave function of the spin singlet 
FQHE has a factor which is precisely this conformal block correlation function. 
It was also natiredea that this factor coincides with the Kalmeyer-Laughlin 
wave function for a Spin Liquid. 

The states with macroscopic spin have a somewhat different physics. The 
existence of a non-zero average field should make a mean field approach more 
sound. The spin sector of this mean field theory has iV T ,x = \N e ±M spin 
up and spin down spinons each feeling an effective uniform magnetic field of 
±t£ M . Because of the SU(2) invariance, there is no Zeeman term and only the 
orbital degrees of freedom see this spin-dependent external field. It is easy to 
see that the highest weight ferromagnetic state with maximal spin is obtained 
by filling up the lowest Landau level of the up spins while leaving the down 
spin sector empty. 

The charge and spin sectors are not decoupled from each other. Firstly, 



the constraint of Eq. 117 sets the local charge density to be the same as the 
local spin density. The wave functions of the allowed states have to satisfy 
this local property. Secondly, if the system is SU(2) invariant, all of the states 
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in a given SU(2) multiplct must have the same filling fraction. Since the 
fully polarized states have to span all of the Jain states for a single-layer 
system, we must conclude that the SU (2) states which are not in a main Jain 
hierarchy cannot achieve the maximum polarization. In other terms, there is 
an upper bound for the spin polarization and, hence, for the total spin itself. 
This argument indicates that the filling fraction v and the spin S of the state 
cannot be completely independent from each other for the allowed states. In 
other words, there should exist a set of selection rules which determines the 
possible combinations of total spin and filling fraction. Similarly it should 
be possible to construct a unified theory of all the FHQE states with SU(2) 
symmetries, instead of the apparently separate descriptions of the spin singlet 
and the fully polarizable states that we discuss here. 

We conclude this Section with a comparison of the states that are obtained 
by this £[/(2)-symmetric approach and the U(l) ® U(l ) theory. A direct in- 
spection of the allowed fractions, Eq. ^ and Eq. 124, for the U(l) ® U(l) 
and SU(2) theories respectively, shows that they do not yield the same al- 
lowed fractions. For instance, the "Fermi Liquid" (compressible) states, with 
the same occupancy of the two layers, allowed by Eq. [79| have filling fractions 
^ = 2,1, §, i, §, |, f , |, §, . . . (r = 1,2,.. .). In contrast, the allowed SU(2) 

"Fermi Liquid" (compressible) states are = 2, |, |, |, |, Clearly, the 

fractions 4 (with k — 1,2,3,...) cannot be realized as SU(2) compressible 
states. Among the states which appear in both hierarchies, we find an in- 
compressible spin unpolarized state at filling fraction jj. It is worth noting 
that there is experimental evidence of an incompressible state at A. It is 
not possible to construct a fully polarized Jain state with this filling fraction 
for single-layers although it may be constructed as a hierarchical state. It is 
strange that this is the only observed fraction for which a hierarchical con- 
struction is needed. It is believed that the experimentally observed state is 
polarized. We have checked that all of the states in levels 1 and 2 in the SU(2) 
hierarchy span the entire level 1 U(l) ® U(l) states. Similarly, the level 4 
SU(2) states span the level 2 U(l) ® [7(1) states, and the levels 3 and 6 SU(2) 
states span the level 3 U(l) <g> U(l) states. However, a large number of in- 
compressible U(l) <g> U(l) states cannot be realized as SU(2) states. In a way, 
this may not be surprising since the U{1) ® U(l) symmetric theory may only 
generate SU (2) as a dynamical symmetry. Nevertheless, it is still surprising 
since in this discussion the form of the interaction terms has not entered and, 
hence, the symmetries of the Hamiltonian have not had a chance to play any 
role yet. However, more surprising is the fact that not all of the SU(2) states 
can be realized as U(l) ® U(l) states. It is easy to check, for instance, that 
the allowed SU(2) state with v = y and S z = has no counterpart in the 
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U(l) ® U(l) states (unless the condition of S z = is relaxed and polarized 
states are also considered). 



5 Ground state wave functions 

In this section we will make the connection between the fermion Chern-Simons 
field theory and the many-body wave-function approach. We will give a de- 
scription in terms of the bilayer system since it already includes the single layer 
system. 

The first step that we need in order to make this connection explicit is 
to find the exact asymptotic form of the density-density correlation functions. 
We need to show first that the long wavelength form of found at this 

semiclassical level, saturates the /-sum rule. This result implies that the non- 
gaussian corrections do not contribute at very small momentum. We will use 
this result to show that the absolute value squared of the ground state wave 
function of this state has the Halperinc3 form at very long distances, in the 
thermodynamic limit. 

The /-sum rule states that the retarded density-density correlation func- 
tion Dqq satisfies 

^ ^-iuD^{u>,Q) = P -QH^ (126) 

This equation implies the conservation of the current in each layer separately. 
It is easy to show that, i n th e basis of the total and relative currents, i. e., in the 
p± = ^ip 2 basis, Eq. |l26| states the conservation of p+ and p_ independently. 
On the other hand, from the general definition of the polarization tensor Kqq , 
it follows that (see for instance reference^]) 

K$(x,y) = -D$(x,y) (127) 



Thus, Eq. 126 also holds if we replace by JTq (except for an extra minus 
sign). 

5.1 Ground state wave function for (m,m,n) states 

We have found that for an (to, to, n) state, the leading order term in Q 2 of the 
zero-zero component of the electromagnetic response is given by 



#00 



a/3 _ P Q 2 ( 1 1 



1 1 



AM to 2 - ((to - n)u>) 2 + ie \ -1 1 



(128) 
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Notice that the poles whose residues in Kqq are proportional to Q 2 , have 
zero momentum frequency given by (m + n)ui — lu c and (m — n)Q. It is 



straightforward to check that Eq. 128 saturates the /-sum rule. Therefore, the 
fermionic Chern-Simons approach gives the correct leading order form for the 
density correlation function, consistent with the /-sum rule, at the semiclassical 
level of the approximation. 

It is important to emphasize that the coefficient of the leading order term of 
Kqq can not be renormalized by higher order terms in the gradient expansion, 
nor in the semiclassical expansion. In the case of the gradient expansion, it 
is clear that higher order terms have higher order powers of Q 2 , and then, 
do not modify the leading order term. In the case of the corrections to Kqq 
originating in higher order terms in the semiclassical expansion, they also come 
with higher order powers of Q 2 . The reason for that is the gauge invariancc 
of the system which implies that higher order correlation functions must be 
transverse in real space or, equivalently, have higher order powers of Q 2 in 
momentum space. Terms in this expansion higher than Q 2 cannot change the 
leading order term. _ ._ . 

We will now follow the methods of references E3E.ll to write the ground 
state wave function in the density representation. We begin by recalling that 
the absolute value squared of the -ground state wave function in the density 
representation | ^ o 1 2 ls given byL3 

l*o[/Oi,p 2 ]| 2 = 

X>A a e _i -f d2x pa< - S) lim (Q\Te l J' d3x A ° ^ ^°^|0) 

(129) 

where Jq(x) = p a (x). 

The operators in this expression are Heisenberg operators of the system 
in the absence of sources. The vacuum expectation value in the integrand of 



Eq. pL 29| can be calculated from the generating functional of density correlation 



functions, Z[A^], defined by 

lim Z [Aq , A a = 0] = / Vi)*V<4)Va a u e «?«-*,iM°,A°) ( 130 ) 

A«(x)^A«(x)<5(x ) J 

Hence, |\l/o[pii P2] | 2 is determined by the generating functional of equal-time 
density correlation functions. 



The path integral on the r. h. s. of Eq. 130 can be written in terms 
of the effective action S e ff(A < ?) for the external electromagnetic field. We 
have seen that, in the thermodynamic limit, and for weak fields, the effective 
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action admits the expansion given by Eq. ^2| Since we need only the density 
correlation functions, it suffice to know the zero-zero component of K?$ . In 
momentum space, and in the small Q 2 limit, Kqq is given by Eq. |128| . Wc 
can see that the dominant term in Kqq is of order 1/B. Higher order terms 
in the gradient expansion will contribute with higher powers of 1/B. The 
same observation applies for all the corrections to Kqq originating in higher 
order terms in the semiclassical expansion. Here the thermodynamic limit is 
crucial since we are only taking into account fluctuations with wavelengths 
short compared with the linear size of the system. The higher order terms, 
which vanish like powers of Q 2 / B, can only be neglected for an infinite system. 

and O, Eq. [l30] becomes 



Using the results of Sections 4.1 



lim Z[A a .A a = 0] = J d2xd2 y A oi s )( lim ^ a ^vo K oo( x >v)) A o 

A%(x)->A%(x)S(xa) 

or, by Fourier transforming the exponent, we get 



"(if) 



(131) 



lim Z[A%,A a =0]=e 



(132) 

The terms dropped in the exponent of Eq. 131 and Eq. 132 represent equal- 



time density correlation functions with more than two densities. These terms 
give rise to three-body corrections terms ( and higher) to the wave function 
that modify the Jastrow form. These higher order terms, as we noted above, 
are of order higher than Q 2 and can be neglected. 

Replacing the expression for Kqq(uj, Q) given by Eq. 128 into Eq. 132 and 
integrating out Aq, we obtain the following form for the absolute value squared 
of the wave function 



\^(xi,... 7 x Nl , yi,...,x N . 



[Xi Xj 



1 2m 



N 2 iVi N 2 

j ii ift-i&rnn 

i<Cj—l i<j—l i=l j—1 

Ni N 2 



n 



B 



x 



exp{ - ^-(^l^l 2 + ^l^l 2 )} 



i=l 



(133) 



where the coordinates Xi are in plane 1, yi are in plane 2, and N% = N2 = — for 
this state. Here we have used that the eigenvalues of the local density operator, 
in a Hilbert space with N a particles, are p a (x) = Yl^i^i® ~ ^*) ~ P°- 
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The wave function of Eq. 133 is the absolute value squared of the Halperin 
wave function c3. Numerical calculations have established E3 that this wave 
function accurately describes the ground state wave function for the (3,3, 1) 
state for d = 1.5£ c , where £ c is the cyclotron length. We have shown that 



Eq. 132 gives the exact form of the ground state wave function at long distances 



and in the thermodynamic limit. Since the leading order term of K q saturates 
the /-sum rule, higher order corrections in the expansion cannot modify this 
result. Arguments similar to the line of reasoning that led to 
used to derive an analogous expression for a single layer systemE 



133 can be 



5.2 Ground state wave function for (m,m,m) states 

For these states, the leading order term in Q 2 of the zero-zero component of 
the electromagnetic response is, according to Eq. ^l| and 94 or ^ 



K 



a/3 
00 



Q 2 



l l 



AM to 2 - lu 2 + it \ 1 1 

p Q 2 

AM w 2 _ y 2Q2 + ie 



1 -1 
1 1 



(134) 



Following the same steps we can prove that the leading order term of Kqq 
saturates the /-sum rule, Eq. 126. All the remarks on the exactness of this 
result are also valid in this case. Substituting Eq. 134 into the expression for 
the absolute value squared of the ground state wave function (Eq. 132) we 
obtain 



|*(fi, ...,x Nl , y u ...,x N2 )\ 2 

Y N 2 

i<j=l i<j=l 

B Ni N 2 



exp{ 



N-l N 2 

»=ij=i 



exp{-^(^| 2 + ^| 2 )} 



Ni 

E 



1=1 

N 2 



Afi N 2 



\X% Xj 



E 2 EE| fi _ f .| I } 



(135) 



where the coordinates Xi are in plane 1, jji are in plane 2, and N\ — N2 = 
for this state. 
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Notice that this ground state wave function is not exactly the same as the 
(to, to, to) Halperin wave function, to which the true ground state approaches 
as d — > 0. There is an extra factor which comes from the fact that there is 
a gapless mode in the spectrum of collective excitations. This contribution is 
analogous to the phonon contribution to the wave function of superfluid He4, 
and just as in that problem, it is essential to obtain the correct properties for 
the spatial correlations of the ground state. This contribution is very small 
at long distances compared to the cyclotron radius, and it can not be written 
only in terms of coordinates in the lowest Landau level. 

6 Conclusions 

In this Chapter, we have reviewed the fermionic Chern-Simons theory of the 
2DEG in a strong magnetic field in the regime of the fractional quantum Hall 
effect. We showed that, for various rational values of the filling fraction v, the 
mean field approximation of this theory describes a set of composite fermions 
filling up an integer number of effective Landau levels, in agreement with Jain's 
picture. This resulting ground state is non degenerate and supports a gap to 
all of its excitations, a property that shares with the true ground state of 
the system. By using standard techniques of perturbation theory one can go 
beyond the mean field approximation in a systematic fashion. 

The theory that we reviewed here is based on a second-quantized fcrmion 
path integral approach. We showed that the problem of interacting electrons 
moving on a plane in the presence of an external magnetic field is equivalent 
to a family of systems of fermions bound to an even number of fluxes, and that 
this theory has the fermions coupled to a Chern-Simons gauge field with Chern- 
Simons coupling constant — (1/27T x 2s). The semiclassical approximation of 
this system has solutions which describe incompressible liquid states, Wigner 
crystals and soliton-like defects. We worked around the liquid-like solution. 
This mean field solution, or Average Field Approximation (AFA), is seen to 
violate Galilean invariance explicitly. At this level of the approximation, the 
center of mass of the system executes a cyclotron-like motion at the effective 
cyclotron frequency. In this sense the gaussian (or semiclassical) fluctuations 
are essential to restore the original symmetries of the problem. We saw that, 
order-by-order in the semiclassical expansion, the response functions obey the 
correct symmetry properties required by Galilean and Gauge invariance, and by 
the incompressibility of the fluid. We showed that, already at the semiclassical 
or gaussian level, the low-momentum limit of the density correlation function 
saturates the /-sum rule, and in that sense this result is exact, i. e., it can not be 
renormalized by higher order corrections. We showed how our methods can be 
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used to calculate the density-current correlation functions and, from them, to 
extract the Hall conductance. We found that it has the correct value already 
at the gaussian level of the approximation. We also obtained the spectrum 
of collective excitations in the low-momentum limit for short-range and for 
Coulomb interparticle pair potentials. 

We have also presented a generalization of our approach to the study of 
the FQHE in double-layer systems and for unpolarized and partially polarized 
single-layer systems. These systems have a richer class of behaviors and excita- 
tion spectra. We have used the fermionic Chern-Simons theory to investigate 
the physics of these states, including results on their response functions and 
for the quantum numbers of their topological excitations. We have also given 
a specific realization of Haldane's criterion for topological stability. 

Finally, we gave an explicit construction of the ground state wave functions 
from our scmiclassical field theory. We show that their universal properties 
are a consequence of general principles, i. e., incomprcssibility and Galilean 
invariance which determine the analytic structure of the equal-time density 
correlation functions at long distances. 
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